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Abstract 



Recent vigorous investigations of topological order have not only discovered new topo- 
logical states of matter but also shed new light to "already known" topological states. One 
established example with topological order is the valence bond solid (VBS) states in quantum 
antiferromagnets. The VBS states are disordered spin liquids with no spontaneous symmetry 
breaking but most typically manifest topological order known as hidden string order on ID 
chain. Interestingly, the VBS models arc based on mathematics analogous to fuzzy geometry. 
We review applications of the mathematics of fuzzy super-geometry in the construction of 
supersymmetric versions of VBS (SVBS) states, and give a pedagogical introduction of SVBS 
models and their properties [arXiv:0809.4885, 1105.3529, 1210.0299]. As concrete examples, 
we present detail analysis of supersymmetric versions of SU{2) and 5*0(5) VBS states, i.e. 
UOSp(N\2) and UOSp{N\A) SVBS states whose mathematics are closely related to fuzzy two- 
and four-superspheres. The SVBS states are physically interpreted as hole-doped VBS states 
with superconducting property that interpolate various VBS states depending on value of a 
hole-doping parameter. The parent Hamiltonians for SVBS states are explicitly constructed, 
and their gapped excitations are derived within the single-mode approximation on ID SVBS 
chains. Prominent features of the SVBS chains arc discussed in detail, such as a generalized 
string order parameter and entanglement spectra. It is realized that the entanglement spectra 
are at least doubly degenerate regardless of the parity of bulk (super)spins. Stability of topo- 
logical phase with supersymmetry is discussed with emphasis on its relation to particular edge 
(super)spin states. 
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1 Introduction 

Strongly correlated systems such as cuprate superconductors, quantum Hall systems, and quantum 
anti-ferromagnets (QAFM) have been offering arenas for unexpected emergent phenomena brought 
about by strong many-body correlation. In particular, the study of QAFM bears the longest 
history since Heisenberg introduced the celebrated quantum-spin model [U [2] and Bethe [3] found 
the first non-trivial exact solution to the quantum many-body problem, and still is providing us 
with attractive topics in modern physics. Generally, in the presence of many-body interaction, it 
is very hard to obtain exact many-body ground-state wave functions and, even if possible, it is rare 
that we are able to write down them in compact and physically meaningful forms. Fortunately, 
in the above mentioned three cases (superconductivity (SC), quantum Hall effects (QHE), and 
QAFM), the paradigmatic ground-state wave functions have been known and greatly contributed 
to our understanding of exotic physics of these systems: the Bardeen-Cooper-Schrieffer (BCS) 
state [4j for SC, the Laughlin wave function |5J for QHE, and for QAFM the valence bond solid 
(VBS) states 0E], which are the exact ground states of a certain class of quantum-spin models 
called the VBS models. In the present paper, we give a pedagogical review of the VBS models 
and their supersymmetric (SUSY) extension, i.e. the supersymmetric valence bond solid (SVBS) 
models [8, 9, 10J, with particular emphasis on their relation to fuzzy geometry. The VBS models 
had been originally introduced by Affleck, Kennedy, Lieb and Tasaki (AKLT) [6] [7] as a class 
of "exactly solvable" models that exhibit the properties conjectured by Haldane [HI [12], namely 
the qualitative difference in excitation gaps and spin correlations in one-dimensional (ID) QAFM 
between half-odd-integer-spin and integer-spin cases. 

As other quantum-disordered paramagnets, the VBS states have a finite excitation gap and 
exponentially decaying spin correlations. However, the VBS states are not "mere" disordered 
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non-magnetic spin states; the spin-S 1 VBS states necessarily have gapless modes at their edges, 
or more specifically, emergent spin-S'/2 edge spins \13\ 114]. This might remind the readers of the 
(chiral) gapless modes at the edge of the QHE systems where excitations are gapful in the bulk 
(see Fig. [I]). Similar features are found in topological insulators as well (EJHSIIIT] and considered 
as a hallmark of the topological state of matter. Furthermore, in ID, the VBS states are known 
to exhibit a non-local order called the hidden string order |18[ 119] . What is prominent in the 
VBS states is that they most typically exhibit a certain kind of topological order even in ID, 
while QHE needs at least 2D to work. This is a great advantage of the investigation of the VBS 
states since in ID most calculations of physical quantities can be carried out exactly by using the 
matrix product state (MPS) representation \20 \ 1211 [22l 1231 124] combined with the transfer-matrix 
method. 

Since the topological character of a system is believed to be encoded in the many-body ground- 
state wave functions rather than in its Hamiltonian, the relatively simple structure of the VBS 
wave function is suitable to investigate its topological properties by using such modern means as 
entanglement entropy [25} [26] or entanglement spectrum [27]. For the above reasons, the VBS 
states, or more generally, the MPS as a class of model states that satisfy the so-called area- 
law constraint [28] have attracted renewed attention as a "theoretical laboratory" in the recent 
study of topological states of matter. In fact, the MPS representation is now regarded as a 
natural and efficient way to describe quantum entangled many-body states, and for a given (ID) 
Hamiltonian the density matrix renormalization group (DMRG) |29| provides a powerful tool to 
find the optimized variational wave function in the form of MPS (see Refs. |30[ \31\ [32] for more 
details about the MPS method and DMRG). The key idea there is that for generic short-range 
interacting systems, only a small part of the entire Hilbert space is important and, depending 
on the problems in question, there are variety of ways to parametrize this physically relevant 
subspace. Although the Wigner's 3j-symbol may give a convenient description of S'C/(2)-invariant 
MPS states [331 EH EH] , we adopt in the present paper the Schwinger-particle formalism (see, for 
instance, chapters 7 and 19 of [36]) to emphasize the analogies to the lowest Landau level physics. 
The list of possible applications includes a convenient description of gapped quantum ground states 
[32 E3 [28] as well as its application to efficient simulations of dynamics [391 SOI El E21 S3 E] and 
variational calculations [IH1 EH H7J 08] . Due to their simple structure, the entanglement entropy of 
the VBS states comes only from the bond(s) cut by the entanglement bipartition [4"9" 1 [50 [ [51 | [52] . 

In some 'anisotropic' MPSs, it is known that a generalized quantum phase transitions (QPT) 
occur as the parameters contained in the MPSs are varied |53 [ I54 [ I55 ] I56]. However, a remark is in 
order about the interpretation of this kind of 'quantum phase transitions' in MPS. Normally, these 
QPTs in MPS are characterized by the divergence of spatial correlation lengths [56J. In generic 
lattice models, on the other hand, the diverging spatial correlation does not necessarily mean 
the vanishing of the excitation gap, while in the Lorentz-invariant systems, these two occur hand 
in hand. In fact, by the structure of MPS, the block (with size L) entanglement entropy never 
exceeds the L-independent value 2 In D [57J (with D being the dimension of the MPS matrix) while 
in quantum critical ground systems whose low-energy physics is well described by conformal field 
theories, the block entanglement entropy is proportional to In L [571 158] . Therefore, one should 
take the QPTs in MPS mentioned above in a wider sense. This kind of phase transitions in the 
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Figure 1: (Color online) Physical analogies between QHE and VBS state. The bulk excitation on 
QHE is gapful while the edge mode is a gapless (chiral) mode. Meanwhile, the bulk excitation on 
the VBS state is gapful while the motion of edge spins is a freely rotating gapless mode. 



VBS-type of states in 2D have been discussed in [33] 133]. 

Finally, we would like to mention the recent application of MPS to the classification of the 
gapped topological phases in ID. As is well-known, there is no true topological phase character- 
ized by long-range entanglement [59]. However, if certain symmetries (e.g. time reversal) are 
imposed, topological phases characterized by short-range entanglement are possible. Since this 
kind of topological phases are stable only in the presence of certain protecting symmetries, they 
are called symmetry-protected topological (SPT) |60j. As has been mentioned above, an appropri- 
ate choice of MPS faithfully represents any given gapped (short-range entangled) ground state. 
Therefore, the problem of classifying all possible gapped topological phases reduces to classifying 
all possible MPSs by using group cohomology (3J3 [B3]. This program has been carried out for 
such elementary symmetries as time-reversal and link-parity in Refs. |6 1 1. [62] (for the fermionic 
systems, see Refs. [64|. 165]) and for the Lie-group symmetries in Ref. [66]. Topological quantum 
phase transitions among these SPT phases have been discussed in, e.g., [60 ] 16? ] 168]. 

In a sense, the second important keyword of this paper, fuzzy (super) geometry obtained by re- 
placing the ordinary (commuting) coordinates with non-commutative ones, is again closely related 
to Heisenberg, who made a pioneering contribution in physics when he had built quantum mechan- 
ics on the basis of non-commuting phase-space variables. Snyder first substantiated Heisenberg's 
idea of non-commutative coordinates in his paper "quantized space-time" |69Q In fact, the VBS 
models have many interesting connections to QHE and fuzzy geometry. To explain the interesting 
relationship among them, let us begin with an analogy between the VBS states and the Laughlin 
wave functions of fractional quantum Hall effects (FQHE). Soon after the proposal of AKLT [6] [7], 
Arovas, Auerbach and Haldane [7T] realized that the Laughlin wave function and the VBS states 
(generalized to higher-spin cases) have analogous mathematical structures upon identifying the 

1 See Ref. [70] for Heisenberg's contribution to the original idea of non-commutative geometry and related historical 
backgrounds. 
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odd integer m, which characterizes the filling fraction v = 1 /m, and the spin quantum number S. 
Specifically, the VBS state is transformed to the Laughlin wave function of the electron system 
on a two-sphere, i.e. the Laughlin-Haldane wave function [72] by using the coherent-state rep- 
resentation of the VBS states and assigning appropriate correspondences between their physical 
quantities. In such translation, the external symmetry (or, S'C/(2)-rotation of the spatial coor- 
dinates) of the Laughlin-Haldane wave function for QHE on a two-sphere is translated into the 
internal symmetry (or, spm-SU(2) symmetry) of the VBS states for the integer-spin chains. This 
analogy can be generalized and one can readily see that the Laughlin-Haldane wave functions 
with some external symmetry are generally transformed to give the VBS states with the identical 
internal symmetry. 

In the past decade, there have been remarkable developments in higher-dimensional general- 
ization of QHE (see Refs. [73J E3j for reviews). So far, the set-up [3 [72] of 2D QHE has been 
generalized to higher- dimensional manifolds, such as 4D [75], 8D [76], S 2n |77] and CP N 
[79] . There also exist q-deformed QHE [80] and QHE on non-compact manifolds [ED E21 E3J 
As is well-known, QHE is a physical realization of the non-commutative geometry [85], and the 
non-commutative geometry brings exotic properties to QHE |86[ \87\ 188] , Therefore, for each 
higher-dimensional QHE, one can think of the underlying higher-dimensional fuzzy geometry, 
such as fuzzy two-sphere [891 EH EI], four-sphere [92], 2n sphere [931 ES EHl ESI E7J [98], CP N 
[991 HOQl EEQ, g-deformed sphere [102] and fuzzy hyperboloids [TO3J QUI [T05l[T06] . 

On the other hand, we have already seen that the Bloch spin-coherent state enables us to relate 
the Laughlin-Haldane wave functions to the VBS states. In correspondence with such QHE, a 
variety of VBS models have been constructed with the symmetries, such as SO(5), SO(2n + 
1) plEl Ell EH], SU{N + 1) P3D1 ED E21 EE], Sp{N) pH EE], and g-deformed SU{2) 
[2D E21 EU EZl EE] [see Fig. |2j. The VBS states demonstrate manifest relevance to fuzzy 
geometry also, by adopting the Schwinger operator formalism (see Sec. 2.1). 

One of the main goals of this review is to illustrate their relations which have been less 
emphasized in previous literatures. As an explicit demonstration of the correspondence, we discuss 
a supersymmetric (SUSY) generalization of the VBS models. Since fuzzy super-spheres have 
already been explored in Refs. pH EI2 EZD E221 E23J and the SUSY QHE in Refs. [EH [1251 
I126[ll27| ^j we can develop a SUSY version of the VBS models (SVBS models) based on the mutual 
relations [51 El [TO]. For the inclusion of fermionic degrees of freedom, the SVBS states exhibit 
particular features not observed in the VBS states with higher dimensional bosonic groups. For 
instance, while the VBS states only exhibit a property of quantum magnets, the SVBS states 
accommodate two distinct sectors, spin sector and charge sector, due to the inclusion of fermionic 
degrees of freedom. Physically, the SVBS states can be interpreted as hole-doped VBS states with 
superconducting property. Mathematically, the SVBS states are regarded as a type of "superfield" 
in terminology of SUSY field theory. As the superfield unifies a various fields as its components, 
the SVBS states realize a variety of ordinary VBS states as the coefficients of the expansion with 
respect to Grassmann quantities. 

By investigating topological features of the SVBS states, we address how SUSY affects the 



2 A variety of super Landau models with super-unitary symmetries have been constructed in Refs. 128, 129, 130 
131, 132, 133, 134 and fuzzy supergeometries have also been investigated in Refs. 135, 1361 11371 1105] . 
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Figure 2: Close relations among fuzzy geometry, QHE, and VBS. They are "transformed" to each 
other with appropriate translations. 

stability of topological phase or entanglement spectrum. As concrete examples, we discuss topo- 
logical feature of the SVBS states with UOSp(l\2) and UOSp{2\2) symmetries, which we call 
type-I and type-II SVBS states, respectively. To perform detail analyses of the SVBS states on 
ID chain, we develop a SUSY version of the MPS formalism (SMPS) [231 EDI EU • MPS formalism 
is now regarded as an appropriate formalism to treat gapped ID quantum systems [371 [28]. Since 
the MPS formalism naturally incorporates edge states, MPS provides a powerful formalism to 
discuss relations between topological order and edge state. Taking this advantage, we can obtain 
general lessons for SUSY effects in topological phases. Specifically, the robustness of the topologi- 
cal phase in the presence of SUSY is discussed in light of modern symmetry-protected topological 
order argument with emphasis on edge superspin picture. We also generalize the lower SUSY 
analyses to the case of higher SUSY; fuzzy four-superspheres and UOSp{l\4) SVBS states. 

The following is the organization of the present paper. In Sec. [2j we give a brief introduction 
to the original SU(2) VBS states and explain its relations to fuzzy two-sphere and the 2D QHE 
through the Schwinger-operator representation and the Hopf map. In Sec.[3j the SUSY extensions 
of fuzzy two-sphere are presented in detail. The corresponding SVBS states and their basic 
properties are reviewed in Sec. |4j In Sec. [5j the SMPS formalism is introduced for the analysis of 
SVBS states. Gapped excitations of SVBS states are derived with use of SMPS formalism within 
single mode approximation. We investigate topological properties of the SVBS states in Sec. [6| 
where the entanglement spectrum and entanglement entropy are derived. (Unpublished results 
for UOSp(2\2) SVBS chain are also included here.) The stability of topological phases in the 



6 



presence of SUSY is discussed too. In Sec. [7j we extend the discussions to the case of fuzzy four- 
superspheres and the UOSp(l\A) SVBS states. Sec. [8] is devoted to summary and discussions. 
In Appendix [Aj we provide mathematical supplements for the super-Lie group, UOSp(N\ 2K). 
Fuzzy four-superspheres with arbitrary number of SUSY are described in Appendix [B| 

2 Fuzzy Geometry and Valence Bond Solid States 

In this section, we give a quick introduction to fuzzy two-spheres and the (bosonic) VBS states. 
Their mutual relations will be discussed, too. 

The original idea of quantization of two-sphere was first introduced by Berezin in '70s |89j . In 
the beginning of '80s, Hoppe [90J explored the algebraic structure and field theories on a quantized 
sphere, and subsequently in the early '90s, Madore [91] further examined their structures coining 
the name "fuzzy sphere". Unlike the ordinary sphere, the fuzzy sphere has a minimum scale of 
area, while respecting the SU(2) continuum (rotational) symmetry as the ordinary two-sphere. 
This is a remarkable property of the fuzzy sphere, when we consider the field theories on it. In 
fact, in the other regularized field theories, such as lattice field theories, the extrinsic cut-off cures 
the UV divergence at the cost of continuous symmetries and the resulting theories only respect the 
discrete space-group symmetry of the lattice on which they are defined. On the other hand, field 
theories defined on the fuzzy manifolds contain an intrinsic "cut-off" coming from the minimum 
area of the fuzzy sphere, and the non-commutative field theories constructed on it were expected to 
have the innate property which might soften the UV divergence to be appropriately renormalized 
in conventional field theories. 

In the middle of '90s, Grosse et al. [92] generalized the notion of the fuzzy two-sphere to con- 
struct four-dimensional fuzzy spheres and supersymmetric fuzzy spheres, i.e. fuzzy superspheres 
\119\ 1120] . In late 90s, the fuzzy geometry was rediscovered in the context of the "second revolu- 
tion of string theory" and attracted a lot of attentions. Researchers began to recognize that the 
fuzzy geometry is closely related to the geometry that the string theory attempts to describe: The 
geometry of multiple D-branes is naturally described by fuzzy geometry (see Refs. [1381 1139] 1140] 
as reviews) and fuzzy manifolds are also known to arise as classical solutions of Matrix theory 
(see e.g. Refs. [93 ] I141j ). Similarly, field theories on fuzzy superspheres provide a set-up for SUSY 
field theories with UV regularization [119} 11421 1143] , and the fuzzy superspheres were also found 
to arise as classical solutions of supermatrix models |144[ 1145] . For details and applications of the 
fuzzy sphere, interested readers may consult Refs. [146} 1147] 1148} [TW\ . Non-commutative geom- 
etry and fuzzy physics also found their applications in gravity [150, 1151} I152j and in condensed 
matter physics [8611871 155]. 

2.1 Fuzzy two-spheres and the lowest Landau level physics 

The fuzzy two-sphere [SU EH [91] is one of the simplest curved fuzzy manifolds [^] The coordinates 
of fuzzy two-sphere, (i = 1,2,3), are regarded as the operators that are constituted of the 

3 The fuzzy spheres also play a crucial role in studies of string theory (see [1531 11541 1155) for reviews.) 
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SU(2) generators 

x, = 2 -fj t . (1) 

Here, R denotes the radius of the sphere and d is the dimension of the SU(2) irreducible repre- 
sentation: 

d = 2j + 1, (2) 

with the SU(2) Casimir index j (J = 0, h, 1, |, 2, • • • ), and Jj (i = 1, 2, 3) are the SU(2) matrices 
of the corresponding representation that satisfy 

[Ji, Jj] = ieijkJk, (3) 

and 

JiJi=jti + l)U, (4) 
where 1<j denotes d x d unit matrix. The 'coordinates' Xi defined in ([!]) satisfy 

[Xi, Xj] = i^e ijk X k , XiXi = R 2 (l- ^\ l d . (5) 

Since J3 takes the eigenvalues j, j — 1, j — 2, • • • , — j, the eigenvalues of X3 are given by 

Each "latitude" of X3 corresponds to a patch of width 2R/d, i.e. the minimum state, on the 
fuzzy sphere^] If the limit of large SU{2) representation d — > 00 (with fixed radius R) is taken, 
the patches become invisible and the discrete nature of the fuzzy sphere is smeared off. In fact, 
one readily sees [Xj,X,] = and XiXi\d-+oa = R 2 and the discrete spectrum ^ of X3 becomes 
continuum raging between —R and R. In this sense, the fuzzy sphere reduces to the ordinary 
(commutative) sphere with radius R in the limit d — > 00. 

As a realization of the fuzzy two-sphere, a convenient way is to utilize the Schwinger operator 
formalism |156|, 1148] and introduce the following operator 

* = (a, bf (7) 

whose components satisfy the following ordinary bosonic commutation relations: 

[a, a f ] = [6, 6 + ] = 1, [a, b] = [a, 6 f ] = 0. (8) 

By sandwiching the Pauli matrices with the Schwinger operators one may simply represent the 
coordinates Xi of the fuzzy two-sphere as: 

X, = ^&ai$, (9) 



4 Here two different viewpoints are possible. One assumes, as is done here, R to be fixed and, in the large-d (or, 
large-j) limit, the width of the patch 2R/d becomes zero to recover the continuous space. The other fixes the width 
2R/d and the radius of the fuzzy two-sphere R ~ d diverges in the large-d limit. The former is reminiscent of the 
scaling limit of lattice field theories upon identifying R with the physical mass. 
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where Ui (i = 1, 2, 3) are the Pauli matrices: 

This is the well-known construction of the 5C7 (2) angular momentum operators introduced by 
Schwinger [157] . In fact, one can easily check that J{ = (\/2)& Oi<& satisfy the standard SU{2) 
commutation relations Q and that the Schwinger-operator representation ^ coincides with the 
original definition ([I]) . The square of the radius is given by 

X l X i = f^j re(re + 2)l n+1 , (11) 

where re denotes the eigenvalues of the number operator for the Schwinger bosons, re = = 
a) a + tfb. Comparing this with eq.((5]), one sees j = re/2 and 

d(n)=n + l. (12) 

Since we utilize the Schwinger operator, the irreducible representation is given by the fully sym- 
metric representation of SU(2): 

\ni, re 2 ) = . = a tWl 6 tW2 |vac), (13) 
Vrei! n 2 ! 

where rei and re 2 are non- negative integers that satisfy n\ + n 2 = re. Physically, |rei,re 2 ) stand for 
a finite number of basis states constituting fuzzy two-sphere which are the eigenstates of X3 with 
the eigenvalues: 

X 3 = ^( ni -n 2 ) = ^(n-2k), (14) 



where k = re 2 = 0, 1, 2, • • • , re. The dimension of the space spanned by the basis states (13) is 



given by d(n). From the expression (14), it is not obvious that the definition ^ respects the 
SU(2) rotational symmetry of the fuzzy two-sphere. However this is an artifact of the particular 
choice of the Xj-diagonal basis states (|13|) . Any other complete sets of the basis states of the fuzzy 



two-sphere can be obtained by applying SU(2) transformations to the basis set (13). In this sense, 
whole Hilbert space of fuzzy two-sphere is "symmetric" with respect to the SU (2) transformation, 
and hence the fuzzy two-sphere is SU{2) (rotationally) symmetric like the ordinary (continuum) 
sphere. 

Another description of fuzzy two-sphere is to utilize the coherent stat^] formalism, which will 
be useful in understanding the connection with the VBS states. The coherent state <f>, or more 
precisely the Hopf spinor, that is labeled by a point on two-sphere (xi, x 2 , X3) {xiXi = 1), is defined 
as a 2-component complex vector satisfying^] 

(xi ■ <Ji)(j)({xi}) = cj>({xi}). (15) 



5 The coherent state here is usually referred to as the Bloch spin coherent state in literature. 
6 The usual coherent 'state' \(j>) (for S = 1/2) defined as \<f>) = &-(f>\Q) satisfies 

{2x v Ji)\ct>) = \4>) . 



Expressing the above in the basis a*\0) and b \0), we obtain eq.(15l 
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In general, the (normalized) coherent state is represented in terms of the Euler angles as 



1 I l + x 3 
y/2(l + x 3 ) Ui + ia?a . 



-5X 



cos ■ 



-5X 



sin 



(16) 



where (xi,X2,x 3 ) = (sin Ocos cp, sin 6 simp, cos 9) (0 < ip < 2n, 0<#<7r, 0<x< 27r), and 
e - 2* denotes an arbitrary C/(l) phase factor. The relation between the point [x\, X2, X3) on a 
two-sphere and the two-component Hopf spinor (it, v) is given by the so-called 1st Hopf mapj^] 



which maps the three-sphere S 3 onto the two-sphere S 2 : 

S 3 A S 2 . 

In fact, a space of normalized two-component complex spinors (ft = (u, v) subject to 



isomorphic to S 3 and eq.(17) gives the mapping onto S 2 

\2 



1. 



(17) 
(18) 

: 1 iS 

(19) 



The Hopf spinor (f> is regarded as the classical counterpart of the Schwinger operator [eq.Q] 
that satisfies 

(20) 



R 



Xi-tfai = -(n + 2)$t. 

If we note that R(n + 2)/d — > R, X$ — > Rxi in the limit d — > 00, the resemblance between ( |15| ) and 
(20) is clear. Multiplying (R/d)& from the right to both sides of Eq.(20) and using 

$t$ 

= n, we 



reproduce Eq.(ll). By replacing the Schwinger operator with the Hopf spinor, one sees that the 



Schwinger operator construction Q of fuzzy coordinates is an operator-generalization (x{ \-¥ Xi 



<p ^ $) of the Hopf map [I55| . 

The fully symmetric SU (2) representation in terms of the Hopf spinor (u, v) is obtained if the 
Schwinger operator in (13) is replaced with the Hopf spinor (a', W) — > (u, v) 



~ I 

yni! n 2 ! 
d_ 

dv ' 



J- 



d_ 
du 



J z 







d 



du dv 



(21) 



with m + ri2 = n, ni,n2 > 0. 

So far, (u, v) are just the two auxiliary variables needed to realize the fuzzy two-sphere. How- 
ever, if we regard them as the physical coordinates of the usual two-dimensional sphere [via the 
Hopf map (18)], we may think that u nin2 represents the wave functions of a certain kind of 



quantum- mechanical systems in two dimensions. In fact, the wave functions u ni m (21) coincide 



with those of the lowest eigenstates of the Landau Hamiltonian on a two-sphere (with radius R) 
in the Dirac monopole background (i.e. the so-called Haldane's sphere [72]): 

1 



H 



2MR* AiAh 



(22) 



7 For construction of higher dimensional Hopf maps, see Ref . |73| for instance. 
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number of \ g g > * J <=c_l_^^> Monopole charge : 

patches: 71 \^=^? ^ II \ S = n/2 

Schwinger-boson S Weyl representation 

construction ^| Wr 

\*J »=^or^ 

Local spin states 

Spin magnitude : S = n/2 



Figure 3: (Color online) One-body-level relationship among fuzzy two-sphere (upper left), Haldane 
sphere (upper right) and local spin states of the VBS state (lower middle). The fuzzy two-sphere 
consists of a finite number of patches, i.e., the basis states, with width 2R/(n + 1). The Haldane 
sphere is a two-sphere with Dirac monopole at its center. The S = n/2 is the monopole charge 
quantized as half-integer or integer by the Dirac quantization condition. In the local spin state of 
the VBS state (lower center), each blob denotes spin- 1/2 degrees of freedom, and n blobs amount 
to S = n/2 local spin by a large Hund coupling on each site. 



where Aj (i = 1,2,3) are the covariant angular momenta 

Aj = -iei jk Xj(d k + iA k ) (23) 
in the presence of the gauge field Ai generated by the Dirac monopole at the origin: 

Ai = 2^mf^y (24) 

In the fuzzy geometry side, the monopole charge n/2 corresponds to the quantized radius of 



the fuzzy two-sphere (11) (see also Fig. |3j). Thus, the lowest Landau level eigenfunctions (on a 
two-sphere) can be "derived" from the fuzzy two-sphere by switching from the Schwinger oper- 
ator formalism to the coherent state formalism (or, the Weyl representation). This is the key 
observation of correspondence between the fuzzy geometry and the lowest Landau level physics 
[73]. 



2.2 Valence bond solid states 

In order to translate the above features to those of QAFM, we first express the spin 1/2 state in 
terms of the Schwinger bosons: 

|t)=a f |vac), ||)=6 t |vac)- (25) 
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The fully symmetric representation (13) constructed out of n Schwinger operators automatically 



realizes the local spin defined on each site with the spin magnitude 



~ ,1 1 In n / \ 

S=(-®-®--.®-) = -. (26) 

Z Z Z miry symm. Z 

Physically, this bosonic construction realizes the ferromagnetic (Hund) coupling among n spin- 
l/2s to yield the maximal spin, S = n/2, at each site. We have already seen that the spin wave 
function, written in terms of the Hopf spinor (u,v), coincides with those of the non-interacting 
electrons moving on a two-sphere in the presence of the monopole magnetic field. 

Now we demonstrate that a similar analogy exists even when we switch on strong interactions. 
Specifically, the exact ground-state wave functions of a class of interacting spin models called 
the valence-bond solid (VBS) model [6J [71] closely resemble the Laughlin wave functions on 
a two-sphere. The ground states of the VBS models (dubbed the VBS states) on a lattice with 
coordination number z (see Fig. [4]) are constructed as follows. As the first step of the construction, 
we prepare n (= z) local S = 1/2 spins {auxiliary spins) on each vertex (site) of the lattice. Next, 
we connect every pair of two spin- 1 /2s on the nearest neighbor sites by spin-singlet bonda^] called 
the valence bonds (VB) (see Fig. [4]). Last, we project the entire 2 n -dimensional Hilbert space of n 
spin-l/2s onto the subspace with the desired value of the {physical) spin S {S < n/2) on each site 
to obtain the spin-5 many-body singlet state. Depending on which irreducible representation we 
use to realize the physical spin, we obtain different states even for the same lattice structure (some 
examples may be found in, e.g., Ref. |109j ). This kind of construction {valene-bond construction) 
applies to any lattice in any dimensions (see Fig. [4] for typical examples in ID and 2D) and 
the state obtained thereby is called the VBS state or, in a modern terminology, the projected 
entangled-pair state (PEPS) [T59l[T60] . 

It is also possible to construct the VBS states for other symmetry groups (e.g., SU{N) [110, 
1111] and SO{N) |107} [108J) with due extension of the notion of the singlet valence bond. In 
these examples, the VBS states thus constructed are constrained by specific symmetry groups 
and normally they do not have any tunable parameteij^J However, if we use the matrix-product 
representation |20j in ID or, in general, tensor-product representation (or, equivalently, the vertex- 
state representation |53t 154] ) in higher dimensions to express the VBS states and consider their 
"anisotropic" extensions of the tensors, it is possible to obtain parameter-dependent states. These 
generalized VBS states are known to have interesting properties. See, for instance, Refs. |20 |I221[55] 
for ID and Refs. [531 15^] for 2D honeycomb- and square lattices. 

Now let us come back to the SU(2) VBS states. If we represent the up and down degrees of 
freedom of auxiliary spin- 1/2 by the Schwinger bosons a) and (see eq.(25)), the singlet valence 
bond on the bond (ij) reads 

{a\b)-b\a)). (27) 



8 Since the spin-singlet state of the two spin-l/2s maximizes the entanglement entropy, the pair is sometimes 
called maximally entangled in modern literatures. 

9 However, the parent Hamiltonians (i.e. the VBS models) may have tunable parameters. For instance, the 
parent Hamiltonian of the spin-2 VBS state (in ID) allows one free parameter to tune even after we fix the overall 
energy scale. 
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Figure 4: (Color online) VBS states on ID and 2D lattices. Filled circles denote auxiliary spin-1/2 
objects which are finally symmetrized to form S = Mz physical spins at each site. Solid lines 
stand for singlet valence bonds between the spin-l/2s. 

Physically, the spin singlet bond denotes the state with no specific spin polarization made of two 
spin 1/2 states, and hence the valence bond represents a non- magnetic spin pairing between the 
two neighboring sites. 

Of course, we can promote the power of the singlet bond from one to an arbitrary integer M 

(4b] - b\a)) M (28) 

to represent M valence bonds between the sites i and j. Thus, one sees that the original valence- 
bond construction [6j [7] of the VBS states is equivalent to the following representation in terms 
of the Schwinger bosons |71j 



|VBS) 



n 

(ij)eN.N 



J 



&tat) M |vac) 



(29) 



where (ij) G N.N. implies that the product is taken over all nearest-neighbor bonds (ij) and |vac) 
denotes the vacuum of the Schwinger bosons. From [Sf + Sj, aj&t — frjaj] = (a = x,y, z), it is 
obvious that the state eq.(29) is spin-singlet. As z bonds emanate from each site of the lattice 
(for instance, z = 2 is for ID chain and z = 2D for D-dimensional hypercubic lattice), we have 

(30) 



Mz Schwinger bosons per site in the VBS state ( 29 ) 



a] 04 + b\bi 



zM, 



and hence the local spin quantum number Si = ^(a\ai + b\bi) is given by 



In particular, for the ID (i.e. 



2) M 



2 1 

Si = -zM. 
2 

1 VBS state, we have 
St = l 



(31) 



(32) 
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Figure 5: (Color online) For the S = 1 VBS state on a finite open chain, there exist spin-1/2 
degrees of freedom at each edge. By construction, the VBS state is the ground state of the VBS 
Hamiltonian regardless of the spin states at the edges. 



and the local Hilbert space is spanned by the following three basis states: 

|1> = ^=af |vac), |0> = 4&}|vac}, |-1) = ^=b?\vac). (33) 



We were a bit sloppy in writing down eq.(29). When considering a finite open chain (with 



length L), we should be careful in dealing with the edges of the system, while, for the VBS states 



on a circle, the expression (29) is correct without any modification. In fact, in (29), the number 
of the Schwinger bosons at the sites 1 and L is M (half of that of the other sites) and we have to 
add the extra edge degrees of freedom represented by the M-th order homogeneous polynomials 
in a) and to recover the correct spin-Ms at the edges. The edge polynomials for M = 1 are 
given by 

/ t (at,6t) = t 9 ^( a t,6t) = 6 t. (34) 

This representation naturally incorporates the physical emergent edge spins |13| [H] localized 
around the two edges. For general M, the precise form of the VBS states on an open chain reads 
as 

L-l 

|VBS) = (a{nb\) M -v J! (a]b] +1 - b)a) +x ) M \a\) u ^$ L )*\s*c) (0<p,q<M). (35) 
Since the VBS Hamiltonian is defined as the projection operator for the Hilbert space of a 



pair of neighboring spins (see Sec. 4.3 for the detail), the state of the form (29) is the ground 
state regardless of the edge polynomials. Therefore, there appear (M + 1) x (M + 1) degenerate 
ground states for the spin-S(= M) VBS state on a finite open chain (see Fig. [5J [7]. From similar 
considerations, it is obvious that the higher- dimensional VBS ground states (e.g., those in Fig. [4]) 
have degeneracy exponentially large in the size of the boundary. In the following, unless otherwise 
stated, we implicitly assume that the VBS state is defined in ID and will focus on the spin-1 case. 

Let us look at some key features of the VBS states in more detail. Classically, AFM on a 
bipartite lattice assumes the Neel-ordered state, where any pair of neighboring spins point the 
opposite directions. To be specific, in the Neel state of (classical) spin-1 AFM, the z-component 
S z of local spins at the site i takes either +1 or —1 depending on to which sublattice the site 
belongs. (Without loss of generality, we may assume that ordered spins are parallel to the z-axis.) 
The spin configuration described by the VBS state is totally different from that of the classical Neel 
state described above. First of all, we note that the VBS state is 5f7(2)-invariant (being a product 
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of singlet bonds between pairs of adjacent sites; see eq.(29) and Fig. ^ and is non-degeneratej^j 
while the classical ground state, i.e. Neel state, is infinitely degenerate with respect to the SU(2) 
rotational symmetry. This implies that though the magnitude of the local spin is Si = 1, its 



expectation value is zero (Sj)vBS = in the bulls and the system is non-magnetic. In this sense 



the VBS state is purely quantum-mechanical and does not have any classical counterpart. 

For better understanding of the VBS state, let us expand the spin-1 VBS state in the S^-basis: 





i 




-i)|vac) 
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1 •••0 + 0- 


•> + l 


••0 + -- 


• • ) + 1 ■ • ■ 00 + • • 


) + 1 • • • 000 


••) 



(36) 

where the coefficient in front of each term on the right-hand side is omitted for simplicity. What 
is remarkable with the above VBS state is that all the states appearing on the right-hand side 
have a very special feature; the states S z = +1 and —1 appear in an alternating manner with 
intervening S z = states. Namely, the ground state exhibits an analogue of the classical Neel 



order called the string order [TH|, [19] albeit "diluted" by randomly inserted zeros (see Sec. 6.1 for 
further detail). Unlike in the case of the classical Neel order, by the SU(2) symmetry of the state, 
the existence of the string order does not rely on the particular choice of the quantization axis 
(z-axis here). 

For the sake of later discussions, we introduce here a concise representation of the VBS state 
and point out remarkable similarity |71| to the Laughlin-Haldane wave function [72] for FQHE 
on a two-sphere. First, we note that by using the SU(2) ~ USp(2) invariant matrix given by the 
2x2 antisymmetric matrix (see Appendix |A[) 



n 2 = i<j 2 =\ , (37) 




the VBS states (29) on generic lattices can be rewritten compactly as 



|VBS)= 11 ($tft 2 $J) M |vac), (38) 
(ij)eN.N. 

where $j denotes the Schwinger operator on the site i 

^ = (44)*. (39) 



10 This is true for a periodic chain and the bulk in a infinite chain. On a finite open chain, the ground state may 
not be spin-singlet and hence may show degeneracy corresponding to the edge states. 

On a finite open chain, (Si)vBS may take non-zero finite value near the two boundaries and decays exponentially 
to zero toward the center of the system. In this sense, magnetism revives near the boundaries and this is the 
manifestation of the emergent edge states. 
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QHE 


QAFM 


Many-body state 


Laughlin-Haldane wave function 

$LH = U^jinVj ~ V iUj ) m 


VBS state 
m=Uh(alb]-b\a]) M \v a c) 


Power 


m: inverse of filling factor 


M: number of VBs between neighboring sites 


Charge 


S = mN/2: monopole charge 


S = Mz/2: local spin magnitude 



Table 1: Correspondences between physical quantities of many-body wavefunctions of QHE and 
QAFM. 



Then, we rewrite the VBS state (38) into the form of the wave function on a two-sphere. 

(40) 



Specifically, by replacing the Schwinger operator with the Hopf spinor 



we obtain the coherent-state (or, Weyl) representation of the VBS state 

\M _ 
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* n 2 



(y>eN.N. 



I! ( u i v j - v i u j) M i 

(ij)eN.N. 



(41) 



Now the formal similarity to the Laughlin-Haldane wave function of QHE 

$lh = H(ti n 2 = Himvj - v iUj r 



(42) 



i<j 



i<j 



on a two-sphere [71] is clearj^] Though the physical interpretation of the quantities appearing in 
these wave functions are different (Table[T]), mathematical similarities between the VBS model and 
QHE may be manifest from the above constructions. Similarities between topological properties 
of VBS and QHE have also been discussed in Refs. [T61~l 1162] . 



2 Precisely, if we use the coherent-state basis 



= . 1 (uja] + Vib]) Mz \v&c) 
for the local spin-Mz/2 states and expand the VBS state in these basis, we obtain the 'wave function' ({ilj}|VBS} oc 

$VBS({«i !«*})• 

13 By the stereographic projection from a two-sphere to a complex plane: 



= v/i 



$lh (42 1 is, in the thermodynamic limit, reduced to the celebrated Laughlin wave function, 

' z i z i 



$ L =n< 



Zi — Zj) e 



(43) 



(44) 
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3 Fuzzy Two-supersphere 



Next, we proceed to the SUSY version of fuzzy two-sphere and generalized VBS states. We 



mostly focus on the cases of the SUSY numbers N = 1 and 2. 
3.1 N=\ 

First, we introduce TV = 1 SUSY algebra, UOSp(l\2) [T551 [Toil fTH5] . which contains the SU(2) 
algebra as its maximal bosonic subalgebra. The UOSp{l\2) algebra consists of the five generators 
three of which are SU(2) (bosonic) generators Li (i = 1,2,3) and the remaining two are SU(2) 
fermionic spinor L a (a = 81,62), which amount to satisfy 

[Li,Lj] = i€ijkL k , [Li,L a ] = -(<Ti)p a Lp, {L a ,Lp} = -(ia2<?i)a.pLi. (45) 

The UOSp(\\2) Casimir is constructed as 

K = LiLi + e a pL a Lp, (46) 

and its eigenvalues are given by j(J + 1/2): j is referred to as the superspin taking non- negative 
integer or half-integer values, j = 0, 1, g, 1, §, • • • ■ The UOSp(l\2) irreducible representation with 
the superspin index j consists of the SU{2) j and j — 1/2 spin representations, and hence the 
dimension of the UOSp{l\2) representation of superspin j is given by 

(2j + 1) + (2j) = 4j + 1. (47) 

The fundamental representation (j = 1/2) matrices of the UOSp{\\2) generators are expressed 
by the following 3x3 matrices: 

h=-h °V i«=U ,° , Tq V (^) 

2 1 ' a 2 \-{ia 2 T a ) o r y ' 



where r a (a = 1, 2) are 



Eq.(48) may be regarded as a SUSY extension of the Pauli matrices. They are "Hermitian" in 
the sense 

l\ = k, 4 = e a plp, (50) 



where J signifies the super-adjoint defined by 

/ 4t r7t\ 

(51) 



C £) / I -St L>t 



'Fuzzy superspheres are also realized as a classical solution of supermatrix models |144l 1145] . 
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Similar to the case of fuzzy sphere ^ , the coordinates of fuzzy two-supersphere [119} I120| are 
constructed by the graded version of the Schwinger operator formalism \V22\ \1'23\ I124j . The graded 
Schwinger operator consists of two bosonic components, a and b, and one fermion component /, 



*= (a, &,/)*, (52) 



with satisfying 



[a, at] = = {/,/*} = !, 

[a, b] = [a, f] = [b, f] = [a, &t] = [a, /t] = [b, /t] = 0. (53) 

21 2 

With use of the coordinates of fuzzy two-supersphere, S* (2 on the left to | denotes the 
bosonic degrees of freedom, while 2 on the right to | does the fermionic degrees of freedom), are 
constructed as 

Xl = -J^k^>, @a = (54) 

which satisfy 

2R R R 

[Xi,Xj] = i—e ijk X k , [Xi,@ a ] = -j-(<Ti)i3a®p, {©q 5 ©/?} = -j(iv2<ri)apXi, (55) 

where d = n + 1 with n = = n-Q + np = a' a + 6^6 + /t/. Square of the radius of fuzzy 

supersphere is given by the UOSp(l\2) Casimir 



( R\ 2 

XiXi + e a p@ a @p = ( — J n(n + i)l 2n +i- 



(56) 



Notice the zero-point energy in (56) reflects the difference between the bosonic and fermionic 
degrees of freedom. The basis states on fuzzy supersphere consist of the graded fully symmetric 
representation specified by the superspin j = n/2: 

\ni,n 2 ) = , = - at ni tf n2 \v&c), (57a) 
Vni! n 2 l 

\mi,m 2 ) = = a) mi b ]m2 /^vac), (57b) 
\Jm\\ m 2 \ 

where nx,n 2 ,mi and m 2 are non- negative integers that satisfy m + n 2 = mi + m,2 + 1 = n. 
|mi,m,2) is the fermionic counterpart of \rii,n 2 ), and thus \ni,n 2 ) and 777-2) exhibit /V = 1 
SUSY. The bosonic and fermionic basis states are the eigenstates of the fermion parity (— l) np 
with the eigenvalues +1 and —1, respectively. Their degrees of freedom are also respectively given 
by 

de = d(n) = n + 1, dp = d(n — 1) = n, (58) 
and the total degrees of freedom is 

d T = d B + d F = 2n + l. (59) 
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The eigenvalues of X$ for (57) read as 



X 3 



(n - k), 



(60) 



where k = 0,1,2,- •• ,2n. Even k (k = 2n2) correspond to the bosonic states (57a), while odd 
k (k = 2wi2 + 1) the fermionic states (57b). Compare the X3 eigenvalues of fuzzy supersphere 
( |60[ ) with those of the fuzzy (bosonic) sphere (14): The degrees of freedom of fuzzy supersphere 
for even k are accounted for by those of fuzzy sphere with radius n, while those for odd k are by 

2| 2 

fuzzy sphere with radius n — 1. In this sense, the fuzzy two-supersphere S { (n) of radius n can 
be regarded as a "compound" of two fuzzy spheres with different radii, and 



71 



Sp 2 (n) ~ S f 2 (i))9S f 2 (n-l). 



1) [Fig.|6j: 

(61) 



The first Sf(n) accounts for the bosonic degrees of freedom (57a), while the second Sf (n — 1) does 
for the fermionic degrees of freedom (57b). 

Similar to the 1st Hopf map, we can construct a graded version of the 1st Hopf map |166|I167] 

V> -> Xi = 2^L^, 9 a = 2^L a ip, (62) 



where U and l a are the UOSp(l\2) matrices of fundamental representation (48), and ip denotes a 
normalized UOSp(l\2) superspinor 



(63) 



with 



ifity = u*u + v*v — 77*77 = 1. 
Here, the super-adjoint of the superspinor is defined by ijr = (u*,v* 



(64) 



-77*) [see also eq.(51)] and 
* represents the pseudo-conjugation. 15 The first two components of tp are Grassmann-even and 



15 The pseudo-conjugation is defined as (77*)* = —r\ and (771772)* = 771772 for Grassmann odd quantities. See 
Ref . [168 for instance. 
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the third component is Grassmann-odd. ip subject to the normalization (64) can be regarded as 

S 2 \ 2 : 

(65) 



a coordinate of the manifold S 3 ' 2 . From (62), we find that Xi and 9 a satisfy the condition of S" 2 ' 2 : 

XiXi + e a p9 a 9p = {ip x ip) 2 



1. 



Consequently, the map (62) represents 



g3\2 _S\^ g2|5 



(66) 



The bosonic part of (66) exactly corresponds to the 1st Hopf map. Note that ip satisfies the 



super-coherent state equation 



Xi ■ kip + e a pl a ip ■ Op = -ip, 



(67) 



and ip is referred to as the super-coherent state or spin-hole coherent state |36J in literature. x% 
are Grassmann even but not usual c-numbers, since the square of x% is not a c-number, xix% = 



1 — CapQaQp, as seen from (65). Instead of X{, one can introduce c-number quantities {yi\ as 

1 



Xi , 



which denote coordinates on S 2 , the "body" of 5 2 ' 2 , as confirmed from 
coordinates on S" 2 ' 2 , ip can be expressed as |124j 



(68) 

1 . With use of the 



iP 



v / 2(l + y 3 )(l+^i^) 



yi + W2 
V(l + 2/3)^1 + (2/1 +iy2)9zj 



( a + *s)(i- 2(rk)^ 
{x 1 + ix 2 )(\ + w ^ ) e 1 e 2 ) | e 

\ (1 + x 3 )9x + (xi + ix 2 )e 2 



(69) 



where e %x stands for the arbitrary U(l) phase factor. The last expression on the right-hand side 
manifests the graded Hopf fibration, S 3 ! 2 ~ 5 2 ' 2 ® S 1 (here, ~ denotes local equivalence): The 
S' 1 (~ [7(l))-fibre, e lx , is canceled in the graded Hopf map (62), and the other quantities, Xi and 



9 a in (69), correspond to the coordinates on S" 2 ' 2 . 



3.2 N = 2 

21 2 

As the geometric structure of S { is determined by the UOSp{\\2) algebra, N = 2 fuzzy super- 

2|4 

sphere S f is formulated by the UOSp{2\2) algebra. The UOSp(2\2) algebra contains the bosonic 
subalgebra, usp(2) © o(2) ~ su(2) © u(l), and is isomorphic to SU (2\1). The dimension is given 
by 

dim[itosp(2|2)] = dim[su(2| 1)] = 4|4 = 8. (70) 
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Denoting the four bosonic generators as Lj (i = 1,2,3) and T and the four fermionic generators 
as L a and L' a (a = 6*1,6*2), we can express the UOSp(2\2) algebra as 

1 



[Li,Lj] — ieij^Lk, [Lj,,L aa ] — -v *) /3a > 
= 0, [r, L aa ] = —e Ta L aT , 



{L acr ,L/3 T } — -S aT (ia2(Ti) a i3Li + -e aT e a /3T, 



(71) 



where L aa = (L a ,L' Q ). Li and L Q form the ?70<Sp(l|2) subalgebra. There are two sets of 
fermionic generators, L a and L' a , which bring M = 2 SUSY. The UOSp(2\2) algebra has two 
Casimirs, quadratic 

JC = LiLi + € a pL a Lp + e a pL a Lp + IT (72) 

and cubic ones [164J. 

To specify a fuzzy manifold, an appropriate choice of irreducible representation is crucial as 
well. The irreducible representation of UOSp(2\2) is classified into two categories; typical repre- 
sentation and atypical representation [1641 1165] . Since the quadratic Casimir (72) is identically 
zero for atypical representation, we adopt typical representation for the construction of 5 f [123J. 
The UOSp{2\2) matrices for typical representation for minimal dimension are represented by the 
following 4x4 matrices: 



1 f <Ti 2 

2 lo 2 2 



V 



2 r Q 

\-(i<W a )* 



2 T a 0' 

-(m 2 T a )* 

0, 

_ 1 /o 2 2 ' 

7 " 2 lo 2 ia 2 



(73) 
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Applying the Schwinger construction, we introduce the coordinates of Sc as 



a 



2R 



where ^ signifies the four-component UOSp{2\2) Schwinger operator 

W = {a, bj, gf, 



(74) 



(75) 



and d = n + 1 with n = = a^a + b^b + y f + g^g. Here a and b are bosonic operators, while 
/ and g are fermionic operators that satisfy 

[a,<J] = [b,tf] = {f,fi} = {g,gi} = l, 
[a, 6] = [a, /] = [a, g] = ■ ■ ■ = {/, 5} = {/, <?t} = . 

2|4 

It is straightforward to evaluate the square of the radius of Sf : 



(76) 



R 



n 2 . 



(77) 
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2|4 

With a given n, the basis states of S* are constituted of the graded fully symmetric representation 
of the UOSp{2\2): 



\ni,n 2 ) ■■ 
|mi,m 2 ) 
\m[,m' 2 ) 

\h,h) = 



\fn\\ n 2 \ 
1 



vac) 



\Jm\\ m 2 \ 
1 



1 



a t mi 6t m2 /t| vac) 
5^ | vac) 



m 



a^&^/^lvac), 



(78a) 
(78b) 
(78c) 
(78d) 



where n\, n 2 , mi, m 2 , m[, m 2 , l\ and l 2 denote non-negative integers that satisfy n\ + n 2 = 
m\ + m 2 + 1 = m'i + m 2 + 1 = l\ + l 2 + 2 = n. The dimension of bosonic basis states, \n\, n 2 ) and 
\h,l 2 ), and that of fermionic basis states, \m\,m 2 ) and \m' 1 ,m 2 ), are found to be equal: 



d B = d{n) + d(n - 2) = 2n, 
dp = 2 x d(n — 1) = 2n, 

and the total degrees of freedom amount to 

c?T = d~B + dp = An. 



(79) 



(80) 



2 1 4 i 

Square of the radius of 5 f (77) does not have the zero-point "energy", since the bosonic and 
fermionic degrees of freedom are canceled exactlj 16 The four sets of the basis states ( 78 ) forms 

^ 2|4 

the M = 2 SUSY multiplet, and suggest the following geometrical structure of 5 f , 



Sf\n) ~ S f 2 (n) Sf(n - 1) S f 2 (n - 1) © Sf(n - 2). 



The latitudes for the the basis states ( 78 ) are given by 

R 



d 



(n - k), 



(81) 



(82) 



with k = 0, 1, 2, • ■ • , 2n. The even ks (k = 2n 2 , k = 2l 2 + 2) correspond to the bosonic basis states, 
(78a) and ( |78b[ ), while odd A:s (k = 2m 2 + 1, k = 2m' 2 + 1) the fermionic basis states, (78c) and 
(78dJ) [Fig. 0. Except for non-degenerate states at the north and south poles X% = ±n, all the 
other eigenvalues of X3 (82) are doubly-degenerate. 



4 Super symmetric Valence Bond Solid States 

In this section, we review the basic properties of the SVBS states [SI |H] and discuss its intriguing 
connection to the SUSY QH wave function and BCS function. 

16 In the cases of fuzzy superspheres based on the UOSp(2\J\f) algebra for J\f > 2, the dimension of fermionic 
basis states is larger than that of the bosonic basis states, and the zero-point "energy" becomes minus. In other 
words, the radius of fuzzy supersphere would become minus for the UOSp(2\N) representations with sufficiently 
small dimensions. Thus, we only deal with the fuzzy superspheres for M = 1, 2. 
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Sf\n) 





Sf(n-1) © 



5 f 2 (n-l) 



— 2l4 

Figure 7: (Color online) Sf is a "compound" made of four fuzzy two-spheres that are considered 
as M = 2 superpartners. The above picture corresponds to n = 2. 



4.1 Construction of SVBS states 
4.1.1 Af = l 

Here, we consider the SVBS states with UOSp(l\2) SUSY (N = 1), which we shall call the type- 
I SVBS states. We apply the mathematical procedure of the constructions of the VBS states 



described in Sec 2.2 The first we prepare is the UOSp(l\2) invariant matrix (see Appendix LAI) 



K 



1|2 



/o 



1 






With the parameter dependent Schwinger operator 



¥(r) = (a,b,Vrf)\ 



the type-I SVBS state is constructed as 



ISVBS-I) 



(ij) 



K 



1[2 



¥j(r)) 



M 



vac) 



UUM 



rflf) 



vac) 



(83) 



(84) 



(85) 



The operators a^, hi and fi are components of the graded Schwinger operator defined on each 
site i and satisfy the commutation relations, [aj,aj] = bjj] = 5{j and {fi, fj} = Sij- Physically, 
the three fundamental states, a^|vac), 6^|vac)and /^|vac), are interpreted as spin f, I and spinless 
hole states (see Table |2 17 Since the fermions always appear in pairs of the form fjft (i, j are 
adjacent), the SVBS states can be regarded as hole-pair doped VBS states, and r stands for a hole 
doping parameter. 

Here, some comments are added. The UOSp(l\2) specific feature does not enter the local 



Hilbert space on each site. For instance, (57) can also be regarded as an irreducible representa- 
tion of 5f7(2|l). Meanwhile, in the construction of the type-I SVBS states (85), the UOSp(l\2) 



invariant matrix was utilized, and then the UOSp(l\2) structure explicitly enters in the many- 
body states. This implies that (super)spin interaction between adjacent sites reduces the SU (2\1) 
symmetry on each site to the lower symmetry UOSp(l\2) in many-body physics. 



In the type-I SVB states (85), the total particle number of the Schwinger particles at each site 
is given by 



zM = a\a i + b\b i + fU i 



(86) 



7 In Ref. [169] . such UOSp(l\2) triplet is dubbed the superqubit. 
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Schwinger operator 


SU(2) quantum number 


Spin state 


at 


1/2 


t) = vac) 


&t 


-1/2 


1} = b* vac) 







\h) = /t|vac) 



Table 2: The physical interpretation of the local states made by the Schwinger operators. /' 
denotes the hole degrees of freedom. 

Since the fermion number f*f takes either or 1, the following two eigenvalues are possible for 
the local spin quantum number Si = \{a\ai + b\bi)\ 

Si = M,M-~. (87) 
In particular, for M = 1, each site can take two spin values 

Si = 1, ~ (88) 
and the local Hilbert space is spanned by the five (4M+1, in general) basis states 

1+2 + ± If 2 

1) = — -a\ |vac), |0) = a\b\ vac), |-1) = — -b] |vac), 

V2 1 11 yjl 1 (89) 

|t) = a\f} | vac), ||) = 6|,/; | vac). 

These constitute M=l SUSY multiplet with the UOSp(l\2) superspin 5 = 1. Similarly, the edge 
states consist of N = 1 SUSY multiplet with the UOSp(l\2) superspin S = 1/2: 

|t)) = at|vac), ||)) = 6t| V ac), |0)) = /t|vac). (90) 

As we will see in Sec(5j the ground state of a finite open chain is nine- fold degenerate (corresponding 
to 3 x 3 matrix-components of the M = 1 type-I SVBS states). 

Intriguingly, the M = 1 type-I SVBS chain interpolates two VBS states in the two extreme 
limits of the hole doping: In the limit r — > 0, |SVBS-I) reproduces the original spin-1 VBS state 
|VBS), 

ISVBS-I) -> |VBS) = Y[(4b\ +1 - £>|a i+ i)|vac), (91) 

i 

while in the limit r — > oo, |SVBS-I) reduces to the Majumdar-Ghosh (MG) dimer state |17U"l 1171] 
|MG), 

|SVBS-I)^|n/l||MG), (92) 



where |MG) is either of the two dimerized states of the MG mode^j [Fig. |8j: 

|MG) = { '* A) = n,even(^! +1 " ^vac), ^ 

i*b) = n i: odd(«l & l+i - bWi+i)\™c 



18 The open boundary condition has been implicitly assumed here; if the periodic boundary condition had been 
used, the two states would have been summed up with a minus sign due to the anti-commutating property of the 
holes. 
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e- 

Figure 8: (Color online) Two MG dimer states related by translation. 

For larger M, |MG) should be replaced with the inhomogeneous VBS states |71] where the number 
of valence bonds alternates from bond to bond. 



In the super-coherent state (63) representatior the SVBS state is expressed as [8j 



*SVBS-I = - ViUj - rrjirjj) M . (94) 

(ij) 



From the Grassmann (odd) property of rj, vI/svbs-i (94) can be rewritten as 



* svbs-i = exp (-Mr V ) • $vbs(K vt}) , (95) 

V ^ UiVj - v iUj J 

where $vbs is the coherent state representation of the original VBS state (41). We can deduce a 
nice physical interpretation of the SVBS states from this expression. Since the SVBS states are 
written as a "product" of the exponential and the original VBS wave function <3?vbs> a h physics 
inherent in the SVBS must be included in this exponential factor: 

expf-MrV ^ )=J]h-Mr ^ ] . (96) 

\ ^— ' UiVj — ViUj I 11 [ UiVn — ViUn I 

Since every time when the factor 

w& (97) 



UiVj — ViUj 

acts to the VBS wave function, the VB between the adjacent sites i and j is replaced with a 
hole-pair (see Fig. [jl: 



UiVj - ViUj > T)iT)j , (98) 



one sees that the SVBS wave function (95) may be expanded as 

* SVBS-I = *VBS - Mr V »± $ V BS + \{Mrf (V W ) 2 + • • • 

^ UiV i+ i - v { u i+ i 2 \^ UiV i+ i - ViU i+ iJ 

+ (-Mr) L / 2 T[ Vi • ( TT " II ) 1 * VB S' (") 

- L .- L 1S - L * UiV i+ l - ViU i+ l 

i «:even i:odd 

Thus, the SVBS chain is expressed as the superposition of many-body states on the right-hand 



side (r.h.s.) of (99). The first term on the r.h.s. is the original VBS chain (which is consistent 



19 In the following discussion, the explicit form of the superspinor is not important. Only the Grassmann property 
of the components matters. 
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with (91)). The second term is the VBS chain with one hole-pair doped, and the third term is 



the VBS chain with two hole-pairs doped. In general, nth term represents the VBS chain with 
(n — 1) hole pairs doped. As the last term, partially dimerized chains, i.e., the VB chains whose 
half sites are occupied with hole-pairs, are realized (see Fig. 10 20 For M = 1, the last term gives 
rise to the Majumdar- Ghosh dimer states, 

(IT IT) ' -•^VBS = (TT - T[)(u i V i+ i-V i U i+1 ) = -$ A + <f> B , (100) 

- LJ - - LJ - U{V j+l — ViUi+l - LJ - - LJ - 

i:cvcn i:odd i:odd i:ovcn 



where $a an d $b are the coherent state representation of |<&a) an d I^b) (93). Now, the physical 
meaning of the SVBS states is transparent: The SVBS states signify a superposed state by all 
possible hole-pair doped VBS states, which can be viewed as a generalization of the resonating 



valence bond state |172j [see Sec. 4.2 for more details]. 



"0—0- ^ -0 



e- 



Figure 9: (Color online) When the exponential factor (97) acts to the VBS state, the factor breaks 
the VB between i and i + 1 site and creates hole-pair instead. (Figure and Caption are taken from 
Ref.0). 

As in the original correspondence between VBS and QHE, the super-coherent state repre- 



sentation of SVBS (94) shows striking analogies to the SUSY Laughlin-Haldane wave function 
[125], 



N 



N 



$slh = n^( r )* Kip ^j(r)) M = JJiiHVj - v iUj - rrHVj) m . (101) 

i<j i<j 

Also for one-particle mechanics, there exist apparent relations between the SVBS and the SUSY 



Landau problem. Indeed, the super-coherent state representation of the basis states (57), 

m( B ) - - „, n U, n 2 



1 

y / rnib"<2! 



(102a) 
(102b) 



gives the lowest Landau level eigenstates of a SUSY Landau Hamiltonian |124j . Here, m, ri2, rn\ 
and ni2 are non-negative integers that satisfy n\ + ri2 = mi + rri2 + 1 = 1. By the stereographic 



projection z\ = Vi/ui and ^ = r/j/uj, the SUSY Laughlin-Haldane wave function $slh (101) is 
transformed to the SUSY Laughlin wave function defined by 



N 



$SL = H(zi - Zj - r^) m e 

i<j 



(103) 



For general M and arbitrary lattice coordination number z, the last term of the expansion realizes a resonating 
valence bond (RVB) state |172] , For instance, on 2D square lattice, M — 2 SVBS state with M = 3 SUSY gives the 
Rokhsar-Kivelson RVB state |173| as the last term. 
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Figure 10: (Color online) The type-I SVBS is a superposed state of hole-pair doped VBS states. 
With finite hole-doping parameter r, all of the hole-pair doped VBS states are superposed to form 
the SVBS state, and the SVBS state exhibits the SC property. At r = 0, the SVBS state is 
reduced to the original VBS state (depicted as the first chain), while r — > oo, the SVBS state is 
reduced to the MG dimer state (depicted as the last two chains). (Figure and Caption are taken 
from Ref.[9]). 



By expanding the polynomial part (zi — Zj — r^j)" 1 in the Grassmann odd quantities, we have 



Z% Zj 



(AT/2)! 



$L, (104) 



where N is the total number of particles (which we assumed to be an even integer) and <J>l 
coincides with the Laughlin wave function on a 2D plane (up to the Grassmann factor e - ^-**^*): 



N 



(105) 



Interestingly, the Pfafnan wave function Pf ( z ^_ z . ) TIi< 7 '( 2; « — Zj) m e ^i ZiZ * for the ground state 



of the non-Abelian QHE [174 J appears in the last term of the expansion (104) at m = 2. 



4.1.2 J\f = 2 

The N = 2 SVBS states, which we call the type-II SVBS states, are constructed as UOSp{2\2) 
invariant VBS states. With the UOSp(2\2) invariant matrix 



n 



2|2 



/ 1 \ 
-10 
0-10 
\ -1/ 



(106) 



and UOSp(2\2) Schwinger operator 



¥(r) = (a, 6, Vrf,^g)\ 



(107) 
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we introduce the type-II SVBS states: 

ISVBS-II) = n(*IW ^2|2 **(r)) M |vae> = JJ( a J&5 - b\a) - rf}f} - rg}g}) M \v aC ). (108) 

(*?} to) 

The inclusion of two species of holes, / and g, allows us to write down a wave function more 
symmetric with respect to the bosonic and fermionic degrees of freedom. The new fermionic 
degrees of freedom gi are interpreted as another species of (spinless) hole, and satisfy the standard 
anti-commutation relations {<?«,<?]} = {fi,gj} = and etc. In the type-II VBS states, there 
appear the local sites such as //<?j|vac) with spin-0, which are not realized in the type-I SVBS 
states. 

We have two species of fermions, and the total particle number of the Schwinger particles at 
each site i reads as 

zM = a\ ai + b\bi + f\fi + g\ gi . (109) 

Since the eigenvalues of nj(i)=fjfi and n g {i)=g\g.i can take either or 1, in the type-II SVBS 
chain (z = 2), the following four eigenvalues are allowed for the local spin quantum number 
Si = l(4ai + b\bi): 

Si = M, M- 1 -, M-±, M-l. (110) 
In particular, for the M = 1 SVBS chain, the local spin values are given by 

^ = 1,^,0, (111) 
and the local Hilbert space is spanned by the following nine basis states, 

1+2 j. x If 2 

I 1 ) = -^f l l vac >> 1°) = a] 6] | vac), |-1) = -j=b\ |vac), 

|t)=al//|vac), ||) = &i/T|vac), 
|f)=aW|vac), ||') = &y|vac), 

|0')=4//|vac). (112) 

The edge states are now given by 

|t)) = a t |vac), |D) = tf\vac), |0» = /^vac), \Cf)) = g*\vac), (113) 

and correspondingly, there appear 4 x 4 = 16 degenerate ground states for the M = 1 type-II 
SVBS chain. 

The M = 1 type-II SVBS chain has the following properties. As in the type-I SVBS state, the 
type-II SVBS state reproduces the original VBS state for r — > 0: 

ISVBS-II) -+ |VBS) = H(a\b\ +1 - 6j a t +1 )|vac) . (114) 

i 

On the other hand, when r — > oo, it reduces to the totally uncorrelated fermionic (F) state filled 
with holes: 

ISVBS-II) |F-VBS) = ±n/k ! |vac). (115) 
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Figure 11: (Color online) Like the type-I SVBS chain, the type-II SVBS chain is also expressed 
as a superposition of the hole-pair doped VBS chains. What is different to the type-I SVBS chain 
is the appearance of the spinless sites, depicted by the large white circles with double holes. The 
MG states are realized in the "middle" of the sequence. The original VBS state and the hole- VBS 
state are respectively realized in the first and last lines. [Figure and Caption are taken from 
Ref.[9].] 



Here we have assumed the open boundary condition Note that, unlike type-I, type-II SVBS 



states have no spin degrees of freedom for r — > oo. The super-coherent state representation of 
ISVBS-II) is given by 

I'SVBS-II = Y[ ( UiV i ~ ViU 3 ~ ~ rr liVj) M 

<ij> 

= exp -Mr } J - -exp -Mr } J -^— ~ -$ V bs- 116) 

By expanding the exponentials in terms of r, the type-II SVBS states can be expressed as a 



superposition of the hole-pair doped VBS states as shown in Fig. 11 



21 If the periodic boundary condition is used, we have zero state for odd-length chains. (The sign factor depends 
on both the parity of the system size and the edge states). 
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4.2 Superconducting Properties 

In both ^svbs-i and ^svbs-II) the fermions always appear in pairs and the wave functions can be 
expressed by a superposition of fermion pairs. We can point out interesting similarities between 
the SVBS states and BCS state of SC fTTSI: 



l BCS > = II h \ J 1 + 9kclcl k )\™c), (117) 
k>0 V 1 + \9k\ 



with electron operator c k and coherence factor g k . |BCS) can be expressed as 
|BCS) oc JJ(l + g fe c[cI fc )|vac) = exp(^# fe 4ct fc )|vac) 



fc>0 fc>0 



|vac) + (^g k c\c ] _ k )\wa,c) + 9kc\c ] _ k ) 2 \va,c) H h g k c\J_ k \vac) . (118) 



2 

fc>0 fc>0 fc>0 



This expansion may remind us of the hole-pair expansion of the SVBS state (99). Furthermore, 
in the limit g k — > 0, the BCS state reduces to the electron vacuum (no fermions), and for g k —> oo, 
it coincides with the filled Fermi sphere | Fermi): 

|BCS) |vac) , (119a) 

|BCS) |Fermi) = Jj4cl fc |vac). (119b) 

fc>0 

These also show apparent similarities with the asymptotic behaviors of the type-II SVBS chain, 
(114) and (115), under the correspondence 

| vac) |VBS), (120a) 

c{cl k o fjgt (120b) 



Doping Parameter r 


' > oo 


Charge-sector 


Insulator Superconductor (SSB) Insulator 


Spin-sector 


Disordered quantum anti-ferromagnets (no SSB) 



Figure 12: The SVBS states exhibit superconducting property in the charge sector with finite r 
in addition to quantum AFM property in the spin sector. 



From the analogies to the BCS state, the SVBS states are expected to exhibit SC property 
in the charge sector by the immersion of hole-pairs to (insulating) VBS states [see Fig. 12 . This 
is quite similar to the mechanism of the Anderson's RVB picture of high Tc SC |172j : A finite 
amount of hole-doping transforms the insulator of resonating valence bond state to high Tc SC 
state. In the following, we explore qualitative arguments for the SC aspect of the SVBS states. 



30 



4.2.1 N = l 

Since the naive SC order parameter such as ( flfl +1 ) vanishes due to the violation of particle 
number conservation at each site, as the SC order parameter of the type-I VBS state we adopt 
the following quantity: 

& = (a k b k+1 -b k a k+1 )fifl +1 , (121) 

which is calculated as 

( A ) = 2 - g • ( 122 ) 

^M{M+l)(l + \r\ 2 ) + \ + \{M+\)) -\{M+\) 2 

It takes the maximum value 



|A_l = (V5-W?^i±^ (123) 



at 



M + 1 )J (124) 
2j]j2M(M + l) y ' 

In particular for M = 1, |r| = | \J — 1-10. The expectation values for the boson and fermion 
numbers, nb(i) = af Oj + &J&i and nf (i) = ///i, are respectively calculated as 

2M + 1 

(n b ) = 2M - 1 + 



V^4M(M + 1)(1 + |r| 2 ) + l' 

2M + 1 , , 

n f = 1 125) 

y / 4M(M + l)(l + |r| 2 ) + l V ; 

As expected, with increase of the hole doping |r|, (nb) monotonically decreases, while (nf) mono- 
tonically increases. The fluctuations for the boson and fermion numbers, Sn 2 ^ = (n^) — (nb) and 
5n 2 = (n 2 ) — (nf) 2 , are given by 

x 2 x 2 2M + 1 (2M + 1) 2 

flref = /)nf = > — ^— — , (126) 

b f v / 4M(M + l)(l + |r| 2 ) + l 4M(M + l)(l + |r| 2 ) + l' v ; 

and their maximum <5nb = <5nf = h is met at 



|r '= 3 ( 1 + 4M(MTT)j- <127> 

The behaviors of such quantities are depicted in Fig. [T3| 

While the BCS state (117) has the particle-hole symmetry 22 the SVBS state does not show 
an exact particle- hole symmetry, r f-> \jr [see the SC order parameter ( 122[ ) for instance]. This 
is because of the unequal properties between boson and fermion operators, (aj&j — bidj) f+ /i/j. 



For (1171, the order parameter Afe = (ctc! fc ), the electron number rife = cj^Cfc, and its fluctuation Srik, are 
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Figure 13: (Color online) Plot of O sc = (A), the hole density (rif) and the hole- number fluctuation 
5nf as a function of r. Here the bulk values are plotted. Inset: Profile of the hole density (r = 0.5) 
for a finite system (L = 20) with different left edge states (t, 4, and "hole"). Only the left edge 
state is changed with the right one fixed to sr =t- The hole density approaches exponentially to 
the bulk value as we move away from the edge. (Figure and Caption are taken from Ref.p].) 



4.2.2 TV = 2 



We define the order parameter for the type-II SVBS state ( 108 ) as follows: 

Aj = (ai&i+a - bia i+1 )(f} fj +1 + g\g\ +l ). 



In Fig. 14, we plotted its expectation value 



<Ai), 



the hole-number (nf), and hole fluctuation dnhoie; 

<n f ) = {fin) = (dm) = K) . 

^ n hole = ( n hole) - (™hole) 2 ("-hole = Uf + Tig) . 



(129) 

(130) 
(131) 



The SC order parameter Osc takes its maximal value at |r| ~ 1.05 for M = 1. Comparing Fig. 13 



and Fig. 14, one may find that the behavior of the order parameter of the type-II SVBS chain is 
more symmetric with respect to \r\ = 1 than that of the type-I SVBS chain. This is because of 
the almost compensation of the contributions from the equal number of the boson and fermion 
species in the type-II SVBS states. 



calculated as 



A fe = 
(nfe) 



t-lk 



i + \g k \ 2 

2 



9k + 91 



\9k\ 



1 + lfffcl 2 ' 
Sn 2 k = ((n k - (n k )) 2 } = 



\9k\~ 
(1 + lfffel 2 ) 2 



9k + g* k ^(fffc+Sfe 1 )' 

Afe and 8nk are symmetric under g k <r¥ 1/gl due to the particle-hole symmetry. 



(128) 
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Figure 14: (Color online) Plot of O sc = (A$), the hole density (whole©) = (fjfi) an d the hole- 
number fluctuation 5n^ ole as a function of r. Inset: Profile of the hole density [r = 0.5] for a finite 
system [L = 20). Only the left edge state is changed with the right one fixed to sr =f- (Figure 
and Caption are taken from Ref.[9].) 



4.3 Parent Hamiltonians 

The VBS state is the exact and unique ground state of a many-body Hamiltonian which we call 
the parent Hamiltonian [HI Ej. The relation between the VBS state and its parent Hamiltonian 
is quite unique. Usually in quantum mechanics, Hamiltonian is firstly given, and then we solve 
the eigenvalue problem of the given Hamiltonian. In most cases, particularly in the presence of 
many-body interaction, it is formidable to exactly solve the eigenvalue problem, and so we need 
to rely on some appropriate approximation method. Interestingly in VBS models, the procedure 
is completely inverse: the many-body state (VBS state) is firstly given, and next the parent 
Hamiltonian is constructed such that its ground state is exactly given by the VBS state. 

We briefly review the procedure for the construction of the parent Hamiltonian. Consider the 
VBS chain with bulk spin 1. The SU(2) decomposition rule of two spin 1 gives the total spin 
J = 0,l,2: 

1® 1 = 00102. (132) 

However, the values for the bond spin of the VBS chain does not take J = 2, since in that case 
all four 1/2 spins are aligned to a same direction and do not form the spin-singlet bond between 



neighboring sites [see Fig. 15 . 



=> 







i ® i J = o e j = i © jX 2 

Figure 15: (Color online) The bond spin J = 2 cannot be realized in the 5 = 1 VBS chain. 
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Hence, the VBS state satisfies the condition 



P J=2 (M + 1)|VBS)=0, (133) 

where Pj denotes a projection operator to the bond spin J = 2. Since the eigenvalue of the 
projection operator is either or 1 for arbitrary adjacent sites, the minimum eigenvalue of the 
"many-body operator", ^ Pj=2(M + 1); is zero. This simple fact is the key observation for the 
construction of the parent Hamiltonian. We then construct the parent Hamiltonian for the VBS 
chain as 

H = ^2Vj =2 Pj=2(i,i + l), (134) 



where Vj=% denotes a positive coefficient 23 It is obvious that the eigenvalues of the parent 
Hamiltonian (134) are semi-positive definite and its ground-state energy is zero. Also remember 
that the VBS state vanishes by the projection operators in the parent Hamiltonian, and hence 

H\VBS) = 0. (135) 



Therefore, the VBS chain is a zero-energy ground state of the parent Hamiltonian (134). Unique- 
ness of the ground state ^ for the parent Hamiltonian (up to the degeneracy coming from the edge 
degrees of freedom) can also be proven [7J. The projection operator is explicitly derived as 

= ,n ./(./.,; ./'(./'.,) • < 136 > 

In the present case 5 = 1, Pj=2(i, i + 1) is given by 

p TT 2S(z)-S(i + l) + 4-J'(J' + l) 
P J=a (.,. + l) = U 6 - J>(J> + 1) 

J'=0,1 v ' 

= \ (s(i) • S(i + 1) + i(S(i) • S(* + l)) 2 + fj • (137) 
After all, for the S = 1 VBS chain the parent Hamiltonian takes the following form: 
H = V £ |s(i) • S(t + 1) + |(S(i) • S(i + l)) 2 + 2 -}, 



(138) 



where the overall proportional factor 1/2 has been absorbed in the coefficient V. The parameter 
V simply determines the energy scale of the system, and is not important in the dynamical 



23 Vj-2 can depend on the lattice site i, but here we postulate lattice translation symmetry and drop the lattice 
index i. 

24 On a finite-size chain, the ground state is, by construction, either unique (for a periodic chain) or degenerate 
with respect to the edge states (for an open chain) . In order to prove the uniqueness in the infinite-size system, one 



has to define the infinite-size ground state carefully [7]- 
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behaviors of the system. The construction of the parent Hamiltonian is rather technical, but the 



resulting parent Hamiltonian ( 138 ) appears to be "physical" : The spin-spin interaction represents 
the Heisenberg AFM, though it contains the quadratic term whose coefficient is one-third of the 
first Heisenberg AFM term. It may be obvious from the above discussion that the VBS state is 



the exact ground state of the parent Hamiltonian, but if the Hamiltonian ( 138 ) was firstly given, 
one might think that it is almost impossible to derive its exact energy spectra and corresponding 
many-body states. However at least for the ground-state, we know the exact form of the wave 
function and its energy. Taking this advantage, we can develop precise discussions about physical 
quantities relevant to the ground-state, such as entanglement spectrum. For excitations, we need to 
rely on some approximation technique to extract useful information from the parent Hamiltonian 



(see Sec. 5.3.2). It should also be emphasized that the present construction can be generalized to 
higher spin VBS states on arbitrary lattice in any dimensions. Generally, the parent Hamiltonian 
is given by 

zM 

h=J2 E vj p j(hi)i ( 139 ) 

(z-l)M<J 

where Vj denotes a positive coefficient and Pj(i,j) is given by 

p T (i i) = TT (s(0 + sQ-)) 2 -^ + i) 
}}j J(J + i)-J'(J' + i) 

_ TT 2S(i)-S(j) + 2S(S + l)-J'(J' + l) 
}}j J(J + 1)-J'(J' + 1) • 1 Uj 

In QHE, the parent Hamiltonian for the Laughlin-Haldane wave function is called the Haldane's 



pseudo-potential Hamiltonian [72]. As discussed in Sec. 2.2 the mathematical structure of the VBS 
state and the QHE is similar to each other, and the pseudo-potential Hamiltonian for QHE can also 
be obtained by applying the translations between QHE and VBS (Table [T]). Indeed, the resultant 



pseudo-potential Hamiltonian for Laughlin-Haldane wave function takes the form similar to (139): 



m(JV-l) 

# = E E vjPj&i)* ( 141 ) 

i<j m(N-2)<J 

where Pj denotes the projection operator to two-body angular momentum J. 
4.3.1 M = l 

By replacing the SU(2) operators with UOSp(l\2) ones, we readily construct the parent Hamilto- 
nian for the type-I SVBS states which are invariant under UOSp{\\2) transformations generated 
by the SU(2) bosonic generators Li, 



L i = «*Hr) [ : J (142) 
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and the parameter-dependent fermionic operators K a 



K a 





(i0"2T a ) 



T, 

1 



a 



*(l/r) 



(143) 



where \P(r) is defined by ( |84| ). Notice that Lj and K a satisfy the f705p(l|2) algebra ( |45| ). With 
use of the C/05p(l|2) generators, it is straightforward to construct the parent Hamiltonian for 
the SVBS states with arbitrary value of the parameter r. First, we need to derive the projection 
operators to the subspaces of bond-superspin J. The decomposition rule for the two superspins 
S is given by 

j = s ® s = o e 1/2 e i e 3/2 e 2 e • • • e 2s. (144) 

Note that the UOSp(l\2) decomposition rule is similar to that of the SU{2) except for the bond- 
superspin decreasing by 1/2. The SVBS state does not contain any UOSp{l\2) bond superspins 
larger than J max = {z — 1)M, and is the exact zero-energy ground state of the parent Hamiltonian 



25 



H, 



type-I 



E E 

(ij) j=j max - 



Vj 



(145) 



where Vj are positive coefficients, and Fj(i,j) are the projection operator onto the superspin-J 
representation of UOSp(l\2) written in terms of the Casimir operator as: 



25 

n 



(K A (t)+K A (j)Y 

J{J+\)-J'(J' + \ 



yr 2K A (i)K A (j) + 2S(S + §)- J'{J' + \) 

J(J+\)-J'{J' + \) 

(146) 

which projects to the two-site subspace of the bond superspin J. Here, K A (i)K A (j) = Li(i)Li(j) + 
e a pK a (i)Kp{j) and K A {i)K A {i) = K A (j)K A (j) = S(S + \). Since the projection operators (146) 
are UOSp(l\2) invariant, the parent Hamiltonian is (145) as well. Following similar discussions 
about the uniqueness of the VBS state in Ref.[7], we can prove that the SVBS state is the unique 
zero-energy ground-state of the parent Hamiltonian (145). 

For concreteness, we demonstrate the derivation of the parent Hamiltonian for the 5 = 1 
SVBS chain (z = 2, M = 1 and then J max = !)■ From (144), we obtain the parent Hamiltonian 



( 145 ) for type-I SVBS chain as 

^{V|P3(i,i + 1) + V 2 F 2 (i,i + l)} 



H, 



chain- 1 



E 



i(5V 2 + 63V,_)K A (i)K A (i + 1) + ^(W 2 - 7Vs)(K A (i)K A (i + l)) 2 
+ ^(V 2 - 5Vi)(K A (i)K A (i + l)) 3 + — {V 2 - 7V3)(K A (i)K A (i + l)) 4 + Vs 



(147) 



Here, we add several comments. Since the Casimir operator (K A (i) + K A (j)) 2 contains pair- 
creation terms of fermions, such as fj fj(&ibj — haj), the Hamiltonian (147) does not preserve the 
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total fermion number Nf = fjfi (though the total particle number, i.e. the sum of boson num- 
ber and fermion number, is conserved). This fermion number non-conserving term is a particular 
structure of the BCS Hamiltonian, and this is in agreement with the SC property of the SVBS 



model (see Sec 4.2). 



The fermionic generators K a are non-Hermitian, and then the type-I parent Hamiltonian ( 145 ) 



is non-Hermitian as well. As a reasonable Hermitian extension of the type-I Hamiltonian, one may 
adopt 

2S 

VjF\{i,j)Fj(i,j). 



H' 



type-I 



£ £ 

(tj) j=j max + 



(148) 



The definition ( 148 ) is a natural generalization of the original type-I parent Hamiltonian, since, if 



the projection operators were Hermitian, from the property 



s2 



Eq.(148) would be reduced 



to the original one ( 145 ) 



4.3.2 M = 2 

The irreducible representation of the UOSp(2\2) group is specified by (j,b), which signify the 
indices of the largest bosonic subalgebra, SU(2) © U(l] 25 The irreducible representation is clas- 
sified into two categories; the typical representation and atypical representation. For b ^ ±j, the 
irreducible representation is called the typical representation with dimension 8j, while for b = ±j 
the irreducible representation becomes to the atypical representation with dimension 4j + 1. The 



eigenvalues of the UOSp(2\2) quadratic Casimir operator (72) for (j, b) representation are gener- 
ally given by 



c = r 



b 2 . 



(149) 



Then for the atypical representation (b = j), the Casimir eigenvalues identically vanish. Since the 
Casimir corresponds to the square of the radius of fuzzy manifold, it may not be probable to ex- 
plore fuzzy geometry based on the atypical representation. Meanwhile, the typical representation 
consists of four SU(2) representations, \b,j,j 3 ), \b + 1/2, j - \, j 3 ), \b- 1/2, j - \, j 3 ), \b,j-l,j s ), 
each of which carries the SU{2) spin index S as 



(\)s = 

(ii) s -- 

(iii) S 

(iv) S 



1/2 
-1/2 
- 1 



(2j + l)-dim , 
2j-dim , 
2j-dim , 
(2j - l)-dim. 



25 We denote the eigenvalues of the SU(2) Casimir operator S 2 as j(j + 1) (j — 0,1/2,1, 
complex number eigenvalues of the U(l) operator B as 6. 



(150) 



.), and arbitrary 
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For instance, 8-dimensional typical representation, (J, b) = (1,0), is constructed with use of the 



components of the UOSp(2\2) Schwinger operator (107): 



1+2 + + If 2 

W 1+) = 2 a i l vac )' 1°) = a l 6 H vac >' R = 2 b i ' vac ^' 



|t) = al//|vac), 



|4.) = 6t/t|vac), 



(151) 



(iii) |t / ) = oM|vac), U') = &W|vac), 

(iv) | ') = « 7 t/t|vac). 

They give rise to = 2 SUSY multiplet. The C/05p(2|2) superspin operators can also be given 
by 

Si = ^(r)Ui^(r), K a = *t(r)i a ¥(r), D a = & (r)Z^(r), B = *t( r ) 7 ^( r ) j (152) 



with li, l a , l' a and 7 (73). Obviously, Si and — iB are the generators of the subalgebra su(2)(Bu(l) 

261 



For a two-site system, the UOSp(2\2) bond superspin operators are constructed as S tot = 
S(i) + SO*). = #a(») + K a (j), Djf = D a {i) + D a (j), B tot = B(i) + B(j), and the quadratic 
Casimir operator is expressed as 



i,j = t > ' 55 + e a/9 A a A /3 ~~ e a/3-L> a Up — tf ti 

= C(i) + C(j) + 2|s(t) • S(j) + e a pK a (i)Kp(j) - e aP D a (i)D p {j) - B(i)B(j) 

= C(i) + C(j) + 2L(i).L(j) 
= 2S 2 + 2L(i)-L(j), 

where C(i) = C(j) = S 2 (S = 0, 1/2, 1, 3/2, • • • ) and L W ' L U) is defined as 

L(i) • L(j) = S(i) • S(i) + e ap K a (i)Kp{j) - e afS D a (i)Dp(j) - B(i)B(j). 

Tensor product of two identical typical representations is decomposed as 

(J,B) = (5,0) ® (5,0) 

— rft25 — 



(154) 



^ 25 - 1 (25 - n, 0) e^o 1 ( 2S ~ 1/2 " 1/2) 



(155) 



(156) 



,25-1 



TOn=0 (25 - 1/2 - n, -1/2) 0f = 7/ (25 - n, 0) . 
For instance pM 

(1, 0) (8) (1, 0) = (2, 0) (3/2, 1/2) (3/2, -1/2) © (1, 0) (1,0) © (1/2, 1/2) © (1/2, -1/2). (157) 



26 The basis states | |l5l| ) carry the 577(2) © Z7(l) indices as 

(i):|6 = 0,j = l>, (ii)±i(iii):|6 = ±l/2,i = l/2), (iv) : |6 = 0, j = 0). 



(153) 



27 The graded fully symmetric representation made of the UOSp(2\2) Schwinger operator carries the SU(2)(BU(1) 
indices, (j,b) = (5,0). 



28 On the right-hand sides of |l57| and (l56j), (1/2,1/2) © (1/2,-1/2) is replaced by not-completely reducible 
atypical representation consisting of semi-direct sum of atypical representations, (0, 0), (1/2, —1/2), (1/2, 1/2), (0, 0) 
(for details, see Sect. 2. 53 in Ref.[168 ). 
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Similar to the U0Sp(l\2) decomposition (144), the UOSp(2\2) bond superspins decrease by 1/2 
[see the right-hand side of (157)]. As also observed in (157), B is specified by J: For integer J, 
B = 0, while for half-integer J, B = 1/2 or —1/2. Hence, the square of B is uniquely determined 
as a function of J: 



B{jf 



1 



(l-(-l) 2J )- 



(158) 



Invoking the usual arguments of constructing the parent Hamiltonian, we derive the type-II 
parent Hamiltonian as 



2S 

-^type-II = ^2 E V J P j( C hj)' 

(l-l/z)2S<J <ij> 

where Pj stand for the projection operators of the bond superspin J: 

Gj - (J' 2 - B{J'f) 



p ^i) = n 



j^j (J 2 - B{J?) - (J' 2 - B(J') 2 ) 



(159) 



(160) 



Here, B(J) 2 and B(J') 2 are given by (158). For L = 1 type-II SVBS chain, the parent Hamiltonian 
is given by 

ffchain-II = ^{^3/2^3/2(^+1) + V 2 P 2 (Ci, i+l )}, (161) 



where the projection operators are 

C M+1 - (J' 2 - B(J'f) 



^3/2(Cm+0 



n 



J'=2, 1,1/2,0 



n 



J'=3/2, 1,1/2,0 



2 - (J 12 - B(J'Y) 

C M+ i - {J' 2 - B(J'f) 
4 - (J' 2 - B(J') 2 ) 



;C 2 i+ i(Ci^ + i - 1)(C- 



4), 



(162) 



96 



C 2 i+i{Ci t i + \ - l)(Ci ji+ i - 2). 



With use of C^i+i = 2L(i) • L(i + 1) + 2, the type-II parent Hamiltonian ( 161 ) can be rewritten as 
^chain-n = EIt^ 2 + 36V V2)(L(i) ■ L(i + 1)) + \(V2 + 3tfy 2 )(L(i) • L(i + l)) 2 

i 

+ ^(5F 2 - 36F 3 / 2 )(£« • ^(i + I)) 3 + \(V2 ~ 12V 3/2 )(L(i) ■ L(i + l)) 4 + V 3/2 \. 

(163) 



It is noticed that, unlike the type-I parent Hamiltonians (145), the type-II parent Hamiltonians 
(159) themselves are Hermitian, since the UOSp{2\2) Casimir itself (154) is given by a Hermitian 
operator. 



5 Supersymmetric Matrix Product State Formalism 

This section reviews the MPS formalism and its supersymmetric version, SMPS. In the formal- 
ism, the edge degrees of freedom are naturally incorporated. Practically, the MPS formalism 
provides a powerful method to calculate physical quantities such as excitation gap, string order 
and entanglement spectrum. 
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5.1 Bosonic Matrix Product State Formalism 



As we have discussed above, the VBS state is expressed as a product of valence bonds defined on 
two adjacent sites. In ID, the VBS state (29) can be rewritten as a product of matrices defined 
on local site [201123]: 

L-l 

|VBS) a/3 = {k 2 $>xt U^\n 2 ^ +1 )(^f L \^) 

= (112^)^1 ■ \J[ K2*t*?J • (^ 2 $^^)|vac) 

= (AiA2---A L ) a p, (164) 
where the "state-valued" matrix A% is given by 



It is clear, by the Schwinger-boson construction, that the row and the column of the matrix A% 
correspond respectively to the valence bonds going from the site i to its adjacent left and right 
sites. Sometimes, it is convenient to write (165) in a slightly different way: 

l 

Ai = A(m)\m)i (representation- (i)) (166a) 



m=— 1 



or 



Ai = A'(a)\a)i (representation- (ii)) 



a=—x,y,z 



\x) = ~ (|+l> - 1-1)) , \y) = ^= (|+l> + l-l)) , \z) = |0) 
In the first representation, the c-number matrices A(m) are given by 



(166b) 



while, in the second, A' (a) is given by the Pauli matrices a a (a = x,y,z). Using these represen- 
tation, we can recast (164) into the form where the c-number coefficients and the basis part are 
separated explicitly: 

{AxA 2 ■ ■ ■ A L ) a/3 = {A(mx)A(m 2 ) ■ ■ ■ A(m L )} a /3 |mi)i®|m2)2<8> • • • ®\mL)L ■ (168) 

{mj} 



The state (164) or (168) represents a collection of the D 2 states (with D being the matrix size) 
specified by the matrix indices (a, f3). In the above case, (a, 0) have a clear physical meaning that 
they specify the states of the two emergent edge spins. 
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It is convenient to represent A{m) [and A* (m)] by the following simple tripod diagrams: 

\A*(m i )\, 3 = .., 1 

a — I j p 



ml 

— □— / 



a, (3 



ml 



(169) 



where the thick- and the thin lines respectively denote the cf-dimensional physical Hilbert space 
(here, d = 3 spin-1 states labelled by m = —1, 0, 1 or a = x, y, z) and the Z)-dimensional auxiliary 
space (D = 2-dimensional space spanned by the spinors a} and , here); the matrix multiplication 
amounts to connecting open thin lines on the adjacent sites. Then, the bra and the ket vectors 



may be depicted by strings of these tripods (Fig. 16 ). Quantum states which can be written in the 

\ m j-i) \mj) \m J+ i) \m L ) 



|MPS) QL ,« R 

= «L 



mi 



m 2 



A A A ,1, A A ... 



(MPS| aL>aR 



3-± j i+1 L 



T I T T 1 1 

(mi| (m 2 | (mj-i\ (mj\ (m L \ 



Figure 16: Diagrammatic representation of MPS and its dual. 



form of ( fl~64l ) or flj~68l ) are in general called matrix-product states (MPS). As has been mentioned 
in Sec. [TJ any gapped gapped short-range states fall into this category. We refer the readers to 
recent readable reviews [31 j. [32] for more details and the applications. 



We would like to comment on an interesting property of MPS (168). When SU{2) rotation 
acts on the state in ( |166b ) as \a) — > Rb a \b), the local MPS matrix Ai transforms like 



At 



SU(2) 



E 

a,b=x,y,z 



R ba A'(a)\b)i= Yl R baO- a \b)i= U j 'a b U\b) i = U^A i U 



(170) 



a,b=x,y,z 



b=x,y,z 



where R is a three-dimensional rotation matrix and U is the corresponding spinor representation. 
Namely, the original SU(2) symmetry for the local spin-1 objects "fractionalizes" into that for 
the two spin-1/2 objects (spinors). From this, it is evident that the spin-1 VBS state on a finite 



open chain [represented by the MPS (164)] transforms, under SU{2) rotation, as if there were two 
spin-1/2 objects ("quark" and "anti-quark") at the ends of the chain. The above is the simplest 
example of more general symmetry fractionalization property of MPS, which will be extensively 
used in Sec. 16.41 

We can generalize the strategy to construct the parent Hamiltonian for the VBS state in 



Sec. 4.3 to any MPS. The idea is to prepare a cluster Hamiltonian and tune the parameters so 
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that the Hamiltonian annihilates all the D 2 states {i.e. matrix elements) of the MPS on that 
cluster. In fact, it can be shown that for any given MPS there exists the parent Hamiltonian 



for which the MPS (164) or (168) gives the (degenerate) ground states |24j . By construction, 
the degree of degeneracy is equal to the number of matrix elements (D 2 ). For example, the 
ground state of the parent Hamiltonian of the spin-S* VBS state (the VBS model) is shown to 
have (5+1) x (5+l)-fold degeneracy, when the model is defined on a finite open chain [3[7j. This 
exact degeneracy on a finite chain is peculiar to the VBS model and, if we slightly deviate from 
the solvable VBS point, the emergent edge spins (S/2) start interacting with each other with a 
coupling constant exponentially small in system size to (partially) resolve the degeneracy. For a 
periodic chain, on the other hand, we have to take the trace over the matrix indices 



|MPS) PBC = Tr j(X)Aj , 



(171) 



and hence the ground state is not degenerate. 
5.2 SMPS Formalism and Edge States 

The Schwinger boson construction described in the previous section can be generalized to SUSY 
cases by using the Schwinger operator which contains both bosons and fermion(s). 

5.2.1 M = l 



Now let us consider the MPS representation for the type-I VBS chain jjjj. The SVBS chain (85) 
is written as a string of 3x3 matrices (a, j3 = 1, 2, 3): 



L-l 



|SVBS-I) Q/3 = (tt 1|2 ¥i(r)*) a J] (■$ i (rf'Jl ll2 y i+1 (r)*)tf*(r)£|vac) 



i=i 



'L-\ 



(^ 1 , 2 * 1 (r)T*iW t - ( Un ll2 %(r)*Mr) f ) • (^ 1 | 2 * L (r)** i (r)t)^|vac) 

(A1A2 ■ ■ ■ Al)al3, 



(172) 



,i=2 



where 



Ai = Ki^i(r)* ■ 1'i(r) t |vac) J 



/ 




(b\) 2 


Mfl 




-(«i) 2 


-4i 


-V~ra\f\ 


V 









/ 


|0>< 


V2\-l)i 


V^\ih \ 




-V2|l>< 


-\0)i 




V 


VF|t>< 







(173) 
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From the expression (172), it is clear that the nine-fold degenerate ground states correspond to 



different possible choices of the edge states: 



isvbs-i: 



open 



Ai 



1=1 



l;s R = t) 

t;*R= t) 
\\sl=°;s r = t) 



\SL= U SR= |) 

\sl= t;«R= I) 
|sl=°;sr= I) 



|«L= I; SR=o; 
\SL= t; SR=0; 

|sl=°; sr=o; 



(174) 



The row index specifies the left edge states and the column one the right. On the left (right) edge, 
the matrix indices {1,2,3} correspond respectively to {|, f, hole} ({f, 4, hole}). 



By looking at the form of Ai (173), one sees that the matrix has a block structure 




F(2,l) 
3(1,1), 



(175) 



where B(m,n) and -F(m, n) respectively denote "bosonic" and "fermionic" (i.e. anti-commuting) 
matrices of the dimension mxn. Therefore, it is convenient to regard Ai as a super-matrix. Thus, 
the SVBS chain can be expressed in the form of supermatrix-product state (SMPS), and the matrix 
size of the SMPS is directly related to the number of edge degrees of freedom. The (S)VBS states 
with different edge states have finite overlaps with each other, which are exponentially decreasing 
as the system size L. That is, two (S)VBS states with different edge states are orthogonal to each 
other only in the infinite-size limit. 

In constructing the SVBS state on a periodic chain, one has to treat the fermion sign carefully 



and one sees that the trace operation used in the standard MPS representation (171) should be 
replaced with the supertrace: 



|SVBS-I) periodic = STY , 

where the supertrace is defined as 

STr(X) = Mu + M22 ~ M 33 . 
From these ^.-matrices, we obtain the following 9x9 T- matrices (transfer matrix): 



A*(aJ)A(a,(3) 





( 1 











2 











r \ 







-1 














































— r 
















-1 



















2 











1 











r 






















r 

























— r 
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\ r 











r 











J 


(a, a, 




= 1 


2,3) 













(176a) 



(176b) 



(177) 
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Here, A* is obtained from A by |-) i— )■ (-| and complex conjugation. 

Using the matrices A(m) and the diagrammatic representation introduced in Sec. 5.1 the 
transfer matrix may be expressed as 



T_. s «=^[^(m)] aJ [Aim) 



m=l 




(178) 



m=l 



The transfer matrix naturally appears in the calculations of the MPS formalism. For instance, 



by using the diagrammatic representation Fig. 16 and the orthogonality of the local basis states 



(m| 

— ^mn; on6 can show that the overlap of the two (S)]VlPSs with, edge states {o^Lo ^r} 
{ckl,ckr} can be written a^ 29 (see Fig. 17 'a)) 

( SLl a R )(MPS|MPS) (QL)aR) = [T L ]^ aL .^ an , (179) 

where the matrix multiplication is taken over the tensor index (a, a). 



ol- □ — □ n-^nn-y-n — n 

(a) ^ ■ ■ ■ m J 

{m l} a L — □ — □ fejsuQyA-^ [fj 



■Q!R 



a,/3 



{ra,} 



□ — □ — □ — □ — □ — □ 



L 

fi — 3 — — 6 — 



R 



Figure 17: (a) Diagrammatic representation of overlap of two (S)MPSs. (b) More general bound- 
ary conditions which are linear combinations of different (ol/r> a L/R,) s may be used. In those 
cases, edge states are expressed by Z) 2 -dimensional vectors 'L' and 'R' (for instance, in the sim- 
plest case (a), the edge-state vector has the components 5a,a L 8 a ,a h )- 



In the periodic case, the above expression is modified: 

(SVBS-I|SVBS-I)pbc = J>gn(a)sgn(/3) {T L } (aM 

a,fi 



where 



sgn(a) 



1 for a = 1,2 
— 1 for a = 3 . 



(180) 



(181) 



For the bosonic MPS, this is straightforward. For the SMPS, one has to treat the fermion sign carefully but, in 
the end of the day, we can check that the final result is the same. 
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Notice that the transfer matrix directly appears in the right-hand side of (180), and hence the 



calculation is boiled down to that of the power of the transfer matrix. The eigenvalues of the 



transfer matrix (177) are computed as 

-l(x3), - ir(x2), *r(x2), ^3 - V%r 2 + 9) , ^ (3 + ^%r 2 + 9) J , (182) 



and plotted in Fig. 
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The largest eigenvalue of the transfer matrix, 1(3 + \/8r 2 + 9 ) , will be 



relevant in the thermodynamic limit. 

10 
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Figure 18: (Color online) Plot of absolute values of the five different eigenvalues of T. The largest 
eigenvalue is always unique and non-degenerate (except for r — > 00). (Figure and Caption are 
taken from Ref.[9j.) 



5.2.2 TV = 2 



Similar to the type-I SVBS chain, we can express the type-II SVBS chain (M = 1) in the form of 
SMPS: 

ISVBS-IIU = (A1A2 • • • A L ) Q p (183) 



where 



At = TZ 2 \2^i(r)^i(ry\v&c)i 



( ottt 

(4? 

rt„t 



-4$ 



Vrfja] -Vrfjb] 



U 



pt„t 



\fra\fl -y/f< 



Mi \ 

r 4a\ 
f\ t 



vac 



/ 



\-VrgJa] -y/rg\b\ rfjgj 
As in the type-I SVBS state, the supertrace is necessary for the periodic system: 

( L \ 



ISVBS-II) = STr 



. 1=1 



(184) 



(185) 
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where STr(M) = M\\ + M22 — -^33 — M44. The corresponding transfer matrix is a 16x16 matrix 
and has seven different eigenvalues: 

{-l(x3), -ir(x4), +ir(x4), - r 2 (x2), r 2 , ^ (r 2 + 3 - /(r)) , i (r 2 + 3 + /(r)) } , (186) 

where /(r) = \A" 4 + 10r 2 + 9. Regardless of the value of r, the largest eigenvalue is ^ (r 2 + 3 + /(r)) 

5.3 Excitations 

In this section, we delve into dynamical properties, i.e., low-lying excitations on the SVBS chain. 



In Sec. 4.3, we have already obtained the parent Hamiltonian from the VBS state. Given the 
form of the Hamiltonian, we can, in principle, investigate dynamical properties of the system. 
Unfortunately, however, even though we have the exact ground state in hand, only limited (exact 
or rigorous) information about the excitations is available |176l 1177] . Nevertheless, when the 
explicit form of the (whether exact or approximate) ground-state wave function is known, the 
single-mode approximation (SMA) gives reasonably good results [87] [TTj : 

VBS state Exact > Parent Hamiltonian SMA > Excitation. (187) 

The SMA not only provides us with a simple transparent way of calculating (approximate) exci- 
tation spectrum but also sets a rigorous upper bound for the true spectrum. 



5.3.1 Fixing parent Hamiltonian 



Since the type-I Hamiltonian (147) contains one extra parameter up to the overall factor, i.e. the 



ratio of V3/2 to V2, we begin with fixing the form of the parent Hamiltonian. One way to fix the 



remaining coupling is to require that the SUSY parent Hamiltonian (147) should reduce to the 



original SU(2) VBS Hamiltonian (138) in the limit r — > 0. This naturally fixes the two coupling 



constants in the type-I parent Hamiltonian ( 147 ) as 



U 3/2 =tanhr, V 2 



V2, 



and we have 



H = V { tanh(r)P t 3 (i, i + 1)P S (i, i + 1) + V^Pj(i, i + 1)F 2 («, * + 1) 



(188) 



(189) 



Some of the matrix elements in the fermionic sector have a factor 1/r, and in the limit r — > they 
are divergent. However, they are harmless in the limit, since the ground states contain no fermion 
in the r — > limit. Consequently, the type-I parent Hamiltonian projected onto the bosonic sector 



coincides with the spin-1 VBS Hamiltonian (138) 



5.3.2 Crackion Excitation 

Now we are ready to derive the excitation spectra by using SMA. The paradigmatic picture of 
the low-lying excitations in the half-odd-integer-spin chains is provided by the so-called Lieb- 
Schultz-Mattis construction |180j . where we apply a slow twist along one of the symmetry axis 
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(say, z-axis) of the spin Hamiltonian. Physically, this boosts the quasi-particles in the system and 
thereby creates low-lying excitations [with energies of the order of ~ (chain length) -1 ] at a special 
momentum determined solely by the total magnetization. Unfortunately, this construction does 
not work in the usual VBS state [23]. Instead, as we will see, an excited triplet bond {crackion 



|181j . i.e., a "crack" created in a "solid" of valence bonds; See Fig. 19) in the VBS gives, to 
good approximation, a physical low-lying excitation. Due to the simple structure of the VBS 
states, the excitations considered in the SMA essentially coincide with the crackions |181] I23j. In 



the Schwinger-boson construction of the VBS states (29), the crackion excitation is obtained by 
replacing one of the singlet valence bonds {a\V- — b\a)-) by a triplet one [either a]at or (ajfrt + b\a)-) 
or &!&!•]• Since a single Schwinger boson o) or &t describes a spin- 1/2 spinon, we may thought of 
the triplet crackion as the confined triplet pair of two spinons. The notion of crackion excitations 
can be generalized in other VBS-type models (states) with higher symmetries [e.g., SU(N)], where 
an intriguing picture on the relation between spinon confinement and the existence of 'Haldane 
gaps' has been proposed [111} 1113] 

Now let us consider the crackions in the SVBS chains [9]. Since we are dealing with a SUSY 



system, we consider two different types of excitations (see Fig. 19) that may be regarded as 
super-partners of each other: 

• Spin excitation: 

Spin triplet excitation (S = 1) created by UOSp{l\2) bosonic operators 

• Spinon-hole excitation: 

Spin doublet excitation (S = 1/2) paired with a hole created by UOSp(l\2) fermionic 
operators 



They are schematically represented in Fig. 19 In the SMA, the (unnormalized) excited-state wave 



function is assumed to be given (for the spin excitation) by 

\k,a) = S a (k)\SYBS) , (190) 

where S a (k) denotes the Fourier transform of the local spin operator S*. Then, the excitation 
spectrum (or, the Bijl-Feynman frequency [178J 1179] ) is obtained by calculating the following 
quantity: 

„ (k,a\n\k,a) (k,a\(H - E )\k,a) 

" smaW = MM" = <MM> ' ( 9 } 



where Ti is given by (189) and E$ is the ground-state energy. We can consider other types of 



excitations by changing S a (k) to other operators. 
5.3.3 Spin Excitation 

Let us start by investigating the action of local spin operators 



S+(i) = a\bi , S-(i) = b\ ai , S z (i) = i(oJoi - b\bi) (192) 
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Figure 19: (Color online) Action of bosonic (spin) operator S z (a) and fermionic generator K\ 
(b) onto the SVBS state. The local operators S a (i) (a = x,y,z) and Ki^{i) respectively create 
a triplet bond and a spinon-hole pair on either of the two adjacent bonds (i — and (i, i + 1). 
[Figure and Caption are taken from Ref.[9].] 



on the SVBS state. A little algebra shows that these spin operators create triplet bonds around 



the site i [see Fig. 19 



S+ISVBS-I 
SflSVBS 



i)=\{-m)+\4 o) )} 



where l^} 1 ^) and \ipl" J ) are obtained by replacing the SUSY valence bond (a\a 



,(o) 



t„t 



t+i 



(193a) 
(193b) 



-rflflx) 



by triplet bonds a|a| +1 and {a\b\ +l + &|a| , 1 ), respectively. By taking the Fourier transform of 



eqs.(193a) and (193b), one immediately sees that the triplon-crackion equivalence (except for the 
momentum-dependent form factor) found in the ordinary VBS states |181|, [23] holds in the SVBS 
case as well. By a simple algebra, it is easy to show the following bound for the true spin-excitation 
spectrum of^ ne {k): 



w sma(^) 



(SVBS-I|S a (fc)(7l - ^ )g Q (-A:)|SVBS-I) 
(SVBS-I|S a (£;)S a (-£;)|SVBS-I) 

1 (SVBS-I| [S a (-k), [H, S a {k)}} [SVBS-I) 

2 (SVBS-I|5 a (fc),S a (-/fc)|SVBS-I) 



(194) 



The last inequality is proven by noting that the left-hand side can be rewritten as the following 
average 

J™duojS aa (k,uj) 



f™duS aa (k, 



(195) 
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and using the spectral decomposition of the dynamical structure factor S aa (k,uj). The spin- 
excitation spectrum obtained [5] in this way is shown in Fig. 20 At r = 0, the dispersion reduces 
to the well-known result of the original VBS chain |71j : 



w sma(^) 



10 
27 



(5 + 3cosfc) (a = x,y,z). 



(196) 



In the limit r /* oo, on the other hand, the spin excitation loses its dispersion. This is easily 
understood by noticing that the ground-state reduces to the Majumdar-Ghosh dimer states on 
which excitations cannot move. 



Figure 20: (Color online) The spin excitation (triplon) spectrum Wsma(^) obtained by SMA. At 
r = 0, it reduces to the well-known dispersion wsma(A;) = 10(5 + 3cos/c)/27 of the spin-1 VBS 
model. When r /* oo (Majumdar-Ghosh limit), dispersion becomes flat. [Figure and Caption are 
taken from Ref. [9].] 



5.3.4 Spinon-Hole Excitations 

The dynamics of doped holes in the spin-gapped background is in its own right interesting [182J. 
In the context of the VBS models, some (both exact and approximate) results have been obtained. 
For instance, in Ref. |183j . the motion of spin-0 holes in the spin-1 VBS background is considered 
and the exact one-hole spectrum is obtained. Motivated by the experiments carried out for hole- 
doped spin-1 compound, Ref. |184| introduced a realistic model and investigated the motion of a 
single spin- 1/2 hole immersed in the gapped spin-1 VBS background. 

A similar strategy can be used to obtain the spectrum of the charged (hole, /') excitation 
which is always paired with the S = 1/2 spinon (a) or 6 T ). These excitations are created by 
applying the two fermionic generators of UOSp(l\2) 



1.0 




K l (i)= l 2 {^_f i a\ + V~rf}b l ) 
K 2 (i) = ±(~f i b\-Vtf t i a i ) 



(197) 
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to the VBS ground state. By using the explicit form of the ground-state wave function, it is easy 
to show 



Ki(i)|SVBS-I) 

,(V2)\ 



l« 2) )-|t (1/2) )} 



(198) 



where the crackion state IVv ) is obtained by replacing the SUSY valence bond {a\b\ +l — b\a\ +1 

•t ' ~ 

i+l 



rf}fj +1 ) with a spinon-hole pair (a\fj +1 + fjol+i) [ see Fig. 19 'b)]. The excited state i^lSVBS) 
is defined similarly with a) in the above expression replaced with ¥. For r = 0, the spectrum is 
given by 



The behavior of the spectrum as a function of r [9] is plotted in Fig. 21 

1.0 




Figure 21: (Color online) The excitation spectrum ui^ MA (k) of a spinon-hole pair obtained by 
SMA. This spinon-hole pair state is created by fermionic generator K\ except at r = 0, where the 
transition matrix elements of K\ from the ground state vanish. [Figure and Caption are taken 
from Ref.[9].] 



We add some comments about distinctions between Wsma(^) an d w sma(^)- Since SUSY relates 
the bosonic generators S and the two fermionic generators K a , one might naively expect the same 
spectra for their corresponding excitations. However, this expectation relies on the existence of a 
'unitary' transformation which linearly transforms the set of the SUSY generators onto themselves. 
Since no such transformation exists in the present SUSY, the spectra for the spin and charge sectors 
indeed exhibit different behaviors. 

6 Topological Order 

As has been discussed in Sec.[TJ no true topological order is possible in ID systems [59]. However, if 
we impose a certain kind of symmetries, there can be topologically non-trivial phases protected by 
the symmetries dubbed symmetry-protected topological phases. One of the typical examples would 
be a non-trivial topological phase with the Majorana edge mode in ID interacting fermions [185J. 
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We already know that there exists an analogous "topological" phase with (almost) free edge spins 
(S = 1/2) at the edges in ID spin systems as well. In this section, we give detailed discussions 
about topological properties of the SVBS states [HI [TO] . In particular, we investigate the string 
order and the entanglement spectrum of type-I and type-II VBS states. Then we generalize the 
MPS argument of the symmetry-protected topological order |61|, [62] to SUSY cases to understand 
the degeneracy structure in the entanglement spectrum. 

6.1 Hidden Antiferromagnetic Order and String Order Parameter 

Before proceeding to the SUSY cases, we briefly recapitulate the hidden non-local order in ID 
spin systems. The concept of hidden order is an isotropic generalization of the Neel order. As 



we have seen in Sec. 2.2, the Neel order for S = 1 antiferromagnetic spin chains looks like (if we 



assume that the ordering occurs in the z-axis) 

••• + - + - + - + -+ •••. (200) 

Here, + stands for S z = +1, and — for S z = —1. Clearly, the spin-spin correlation 

(S?S*) (201) 

depends on the parity of the number n of sites between i and j. A simple way of turning this 
position-dependent (or, alternating) correlation into the smooth ferromagnetic one would be to 
insert a phase factor (— l) n between the two spins 

(5|5| + J i y {S z (-l) n S z +n ) (202) 

in order to cancel the sign factor coming from the alternating +1 and —1 between the sites j and 
j + n. Namely, the Neel antiferromagnetic correlation translates to the ferromagnetic correlation 
in (5|(-l)«5| +n ). 



On the other hand, as is seen in the expansion of the VBS state [see (36 )], a typical S z sequence 
appearing in the state reads as 

••• + - + 0- + - 00 + - 0+ (203) 



As has been pointed out already in Sec. 2.2 by removing zeros in the sequence, we can reproduce 
the usual Neel order. In this sense, there still exists a certain kind of Neel order, though "dis- 
ordered" by randomly inserted 0s, called the hidden string ordeip^j \18\ I19j . However, the trick 
used above does not work since, due to the intervening 0s, the positions of +1 and —1 are random 
(though they still appear in an alternating way) and the phase (—1)™ can not cancel the sign 
factor. Nevertheless, a little thought tells that the following choice will do the job: 

exp {mS z tot {j, j + n)) = {-!)* of ±x between j and j + n , (204) 



30 However, there is a striking difference from the usual Neel order. As has been mentioned in Sec. 2.2 the string 
order exists regardless of the choice of the quantization axis, while the Neel AF order is observed only in a particular 
direction. 
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where 

SUhj) = J2 S k ( 205 ) 

k=i 

stands for the partial sum of Sf. between an arbitrary pair of sites i and j. Therefore, it is 



suggested that we should use, instead of eq.(202), the following pair of non-local order parameters 



^t'ring = lim ( S* exp{mS* t (j + 1, j + n)}S* +n 

j I (206) 

«fflg= I™ ( S* exp{mSf ot (j,j + n-l)}S| +n ) 



known as the string order parameters |18(, [19] in order to characterize the non-trivial spin order in 
the (S = 1) VBS state. In the first line, S x ot (j + 1, j + n) is defined similarly to S* ot (j,j + n — 1) 
[see eq.pOo]]. 



The above expressions have been guided by a simple physical intuition. However, as has been 
pointed out in Refs. |186[ I187j . the string order parameters have in fact a deeper meaning than 
we expect from the above simple argument. To see this, first we note that the string correlation 
functions can be recasted in a suggestive way: 



Ol^Ahj + n) = ( Sf exp{i7r5f ot (j + 1, j + n)}S?_ 



j+n 



QX OX 



OrtxingUd + n) = {S z j exp{ivr5 t 2 ot (j, j + n - l)}S z j+n 



(207) 



where S x and S z are defined as 

S* = S% exp{mS? ot (j + 1, L)} , SJ = exp{^ ot (l, j - 1)}S Z . (208) 

The physical meaning of these operators may be best understood by considering the 'height plot' 
of the VBS spin configurations [TBI [23], where the value S z ot (l,j) is represented by the "height" 



between the sites j and j+1 (hence 5| itself is a "step" at the site j). Fig. 23 'a) is a plot of a typical 
S^-configuration of the spin-1 VBS state when the left edge state is y. One can clearly see that the 
meandering steps are always confined between the heights and +1. A similar analysis in the case 
with the left edge state \, shows that the heights are either or —1. From these observations, it 
is concluded that the string exp{z7rSf ot (l, j — 1)} attached to the left of the operator SJ somehow 
suppresses the strong fluctuations in S z and that S z takes either or +1 (0 or —1) when the 
left edge state is f {\). Similarly, one can show that S x becomes weakly ferromagnetic depending 
on the right edge states. In short, non-zero string order parameters translate to the existence 
of a certain kind of weakly ferromagnetic order in the x and the z directions (i.e. (S a ) ^ for 
a = x,z). 
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Remarkably, the following unitary transformation |188j 

U KT = exp i in ^ S z k S* I = ]J exp{inS^} , (209) 

[ k<3 J k<j 

relating the two operators S a and S a as 

Sj = UkT Sj ^kt ' = ^ KT ^KT i (210) 

transforms the original [SU(2)-invariant] Hamiltonian into the one ?7kt H ^kt which is invariant 
only under the dihedral group D2 (or Z2XZ2) consisting of two tt rotations with respect to the 
x and the z axes |186|, I187j . Therefore, we can interpret the existence of the string order (both 
in x and z) in the original system as a consequence of the spontaneous breakdown of the Z2XZ2- 
symmetry in the transformed system Ukt'H ^kt anc ^ * ne resulting (weak) ferromagnetic order in 
S x > z . 

The idea of non-local hidden order and edge states has been to some extent generalized [162, 
1188} 1231 1189] to other values of integer-spin- S although the hidden Z2xZ2-symmetry is never 
broken [188] in the case of even-S (see Fig. 22)^J In the course of these studies, it has been 



recognized that there are some differences [188, 231 in the ground-state properties according to 
the parity of S. Nevertheless, by analogy with the quantum-Hall systems |161j . the ground state of 
generic integer-spin antiferromagnetic chains, including the original VBS state and its higher-spin 
generalizations [71], characterized by certain kinds of non-local correlations and emergent edge 
states have been called 'topological' in a rough sense. 

To illustrate distinct behaviors depending on the parity of spin S, we may introduce the 
following generalized (angle-dependent) string order parameter [188J: 

CffigW = }%( S j exp{^ ot (j, j + n - l)}S* +n \ , (211) 



where the parameter 9 has been introduced for convenience When 9 = 0, it reduces to the 



usual spin-spin correlation and, when 9 = ir, it coincides with the string order parameter discussed 



above. The behaviors of the generalized string order parameter (211 ) are shown for several values 



of bulk spin S in Fig. 22 [23]. The string order parameters are symmetric with respect to 9 = tt 
and generally have S peaks for the spin-S 1 VBS state. As is expected from that the ground state 
is magnetically disordered, the infinite-distance limit of the usual spin-spin correlation function 



0) vanishes regardless of the value of S. On the other hand, Fig. 22 demonstrates that 



^string takes finite values for the odd-S 1 VBS states while it vanishes for even-S. Therefore, in 
the sense of the Z2xZ2-symmetry argument [1861 H87] mentioned above, this hidden symmetry is 
never broken |188| in the even-S" VBS models. 



31 This does not mean that Z2 xZ2-symmetry never breaks down in any even-S chains. In fact, even when S =even, 
one can construct the model ground states which have non- vanishing string order parameters. 

32 The introduction of the 6 parameter is mainly motivated by the idea that the intermediate string might somehow 
cancel the fluctuations between the two spins S z (i) and S z (i + n). For 6 = 7T, this works perfectly in the S — 1 VBS 
state. Except for 9 = (ordinary spin-spin correlation) and 6 — tt (string correlation), no symmetry-related reason 
has been found so far. 
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6/tt 

Figure 22: (Color online) The behaviors of string order parameters in correspondence with magni- 
tude of bulk spins. In particular, at 6 = ir, the string order parameters of even spin VBS (S = 2, 4) 
vanish, while those of odd spin VBS (S = 1,3) have finite values. 



6.2 Generalized Hidden String Order in SVBS Chain 
6.2.1 Af = l 

Unlike the original VBS chains, S z = 1/2 and —1/2 generally appear in the S z sequence of the 



SVBS chain, and a typical S z sequence of the SVBS chains is given b} 33 

• • • t t II +-00 t I + £J £^ • • • . (212) 

From the sequence, one can "derive" the ordinary hidden oder. First, we search the spin-half sites 
from the left, and whenever we encounter a pair of spin-half sites we sum the two S z values to 
have the effective S z (= +, 0, — ) (e.g. | |i— )• — ): 

•••0 + 00 - + - 000 + 0- 0- -- . (213) 

Next we remove the zeros in the sequence to recover the standard Neel pattern: 

••• + - + - + -•••. (214) 



This observation leads to the existence of (generalized) hidden order of the SVBS chain. Fig. 23 
shows a typical height configuration corresponding to the usual VBS chain and its SUSY coun- 
terpart. Reflecting the existence of hidden order, the height configuration is always meandering 



between the height and the height 1 34 It should also be noted that the height-configuration 

33 It may be worthwhile to give some comments on the relation to the ferrimagnetic spin chains that also consist 
of alternating spin 1 and spin 1/2 |46| . Though both the ferrimagnetic chains and the present SVBS chains contain 
spin-1 and spin- 1/2 degrees of freedom, in the SVBS chains the spin 1 and spin 1/2 are not necessarily alternating 
[see ( |212[ ) and Fig. [To] . More importantly, while ferrimagnetic spin chain can exhibit a long-range magnetic order 
as the order parameter commutes with the Hamiltonian, the ground state of the SVBS chain itself is spin-singlet 
and the S(7(2)-symmetry is never broken (spontaneously). 

34 The same reasoning applies to the general spin-S 1 VBS cases and the height configurations are confined within 
a region of width S |23| . 
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Figure 23: (Color online) Height plot of typical spin configurations in S = 1 VBS chain (a) and 
5 = 1 SVBS chain (b). Note that heights are confined within a region of width 1. Although a 
simple "diluted" Neel picture does not hold because of the presence of hole pairs, still we can find 
string order when hole pairs are grouped together in (b). (Figure and Caption are taken from 
Ref.[9].) 



(205) is directly reflected in the components of the SMPS (174) 
/ 
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To substantiate the existence of the hidden order, we explicitly calculate the string parameter for 



the SVBS chains. The behaviors of the string order are depicted in Fig. 24 with respect to the 
hole doping parameter. 
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Figure 24: (Color online) The string correlation function 0^ ring (206) of the SVBS infinite-chain 
for several values of the superspin 5 = M is plotted as a function of r. Notice that, in the limit 
r — > 0, the string order parameter 0^ ring for the S = M SVBS chain reproduces that of the 
5 = M VBS chain. (Figure and Caption are taken from Ref.[9].) 



From Fig. 24, one can find distinct behaviors of the string parameter in terms of the parity 
of bulk superspin. The string parameters of odd superspin SVBS chains generally decrease with 
increase of the hole doping, while those of the even superspin SVBS chains increase. Since the 
hole-doping simply reduces the spin degrees of freedom on the spin chain, the decrease of the 
string order of the odd superspin SVBS chains may be naturally understood. On the other hand, 
the string order behavior of the even superspin SVBS chains is quite interesting, since the string 
order revives with the hole doping. Intuitive explanation may go as follows. For instance, consider 
5 = 2 SVBS chain. At r = 0, the 5 = 2 SVBS chain is exactly identical to the S = 2 VBS chain 
which is essentially constituted of two 5 = 1 VBS chains, upper and lower chains. By doping the 
holes to the 5 = 2 VBS chain, the spin degrees of freedom, say, on the upper 5 = 1 VBS chain 
decrease to have deficits in the chain. Below the deficits on the upper chain, the spin degrees 
of freedom on the lower 5 = 1 VBS chain emerge and come into effect. In this way, the spin 
degrees of freedom of the lower 5 = 1 VBS chain contribute to generate the finite string order 
with increase of hole doping to the 5 = 2 VBS chain. This intuitive explanation can be applicable 
to the revival of the string orders of general even superspin SVBS chains. Consequently, the SVBS 
states bear finite string order with a finite amount of hole-doping regardless of the parity of bulk 
superspin. This is the salient SUSY effect to the topological stability of quantum spin chains. We 



revisit this effect in the context of the symmetry protected topological order in Sec 6.3 and 6.4 



6.2.2 M = 2 

For the 5 = 1 type-II SVBS infinite chain, the string correlation is computed as 

oo 4 
^Tring = ( r 2 + 1) p + 9 ) ' ( 216 ) 
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which is plotted in Fig. 25 The string order (216) takes the value of the S 
at r = like the S 



1 VBS chain, 4/9, 

1 type-I chain, while approaches to zero in the r /* oo limit unlike the 5 = 1 
type-I chain. Since spin degrees of freedom completely disappear from the type-II chain in the 
r / oo limit [Fig. [TT] , the string order of the type-I and II chains shows qualitatively different 
behaviors. 




1 + r s 



Figure 25: The infinite-distance limit of the string correlation function ( 206[ ) as a function of r. 
The value of string correlation smoothly decreases from the original VBS value 4/9 to (no spins 
left). (Figure and Caption are taken from Ref. |9j.) 



6.3 Entanglement Spectrum and Edge States 

As discussed above, the string order of type-I SVBS chains with even bulk-superspin S revives upon 
hole doping. This suggests that, in contrast to their bosonic counterpart, even the SVBS states 
with even superspin can host the stable topological phase. Though the string order parameter is 
appealing in its similarity to the order parameter in FQHE |161j and its relation to the hidden 
Z2 xZ2-symmetry, its fragility under perturbations has also been discussed recently [190} 160] and 
the alternative "order parameter" has been sought for. 

Recent development in quantum-information-theoretic approaches to quantum many-body 
problems enables us to extract information on the bulk topological order from the entangle- 
ment properties of the ground-state wave function [19H 11921 127j . The topological states in one- 
dimensional (ID) spin systems have been reconsidered [60"1 16 1 (. [62] from the modern point of view 
and the precise meaning of the topological Haldane phase has been clarified. In these studies, the 
string order parameters and the edge states, which in general are not robust against small per- 
turbations, are replaced by more robust objects (i.e. the structure of the entanglement spectrum 
or the structure of tensor- network). In particular, it has been shown in Ref. |61[ W2\ that the 
existence of (at least one of) the discrete symmetries (time-reversal, link-inversion and Z2 x Z2 
symmetry) divides all states of matter in ID into two categories: topologically-non-trivial ones 
and the rest. Generic odd- integer-,!? spin chains belong to the former while even-S chains to 
the latter. The hallmark of the topological phase protected by the above discrete symmetries is 
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that all entanglement levels are even-fold degenerate. In this formulation, the difference between 
the odd-S* VBS states and the even-5 ones is naturally understood in terms of the entanglement 
structure; the degenerate structure exists only for odd-S caset 35 It should also be mentioned 



that the topological phases of one-dimensional gapped spin systems have been classified by group 
cohomology, [591 E3] an d the detailed analyses based on the Lie group symmetries are reported in 
Ref. (66]. 

Following the proposal of Li and Haldan^j |27j . we use the structure of the entanglement 
spectrum (e.g. degeneracy of the entanglement levels) as the fingerprint of topological phases. 
Then, the problem of the topological stability of the SVBS chains translates to the stability of the 
degenerate structure of the entanglement spectrum. Before proceeding to the details, we briefly 
introduce characteristic features of the entanglement entropy of the original S'C/(2)-invariant VBS 
states. It has been reported that the entanglement entropy of the SU (2) spin-S" VBS state on an 



infinite chain is given by a constant 37 determined essentially by the degrees of freedom of the edge 
spins 5/2 [501 EI]: 

Se.e. ~ log(5 + 1) (L -> oo) . (217) 

Corresponding to the parity of the bulk spin of the VBS chains, there appear either integer or 
half-integer spin at the edge. Meanwhile, in the presence of SUSY there necessarily appear both 
integer and half- integer spins at the edge, since SUSY relates integer and half-integer spin degrees 



of freedom [Fig. 26 . Such particular feature of the edge spins is crucial in understanding of the 
salient structures of the entanglement spectrum of the SVBS chains. 



Bulk spin 



Integer - 



Edge spin 

Even 4 ► Integer - 

SUSY | 

Odd 4 ► Half-odd-integer - 



■> Half-odd-integer 

Figure 26: (Color online) SUSY relates the edge spins with different parity. 



6.3.1 Schmidt Decomposition and Canonical Form of MPS 

Before going into the detailed discussion, it would also be worthwhile here to give the derivation 
of the entanglement spectrum using MPS. Suppose we divide a system into the two parts A and 

35 This does not mean that all the odd-S 1 spin chains have the degenerate entanglement spectrum. We can construct 
an odd-S 1 spin state without the even-fold degeneracy. 

36 Li and Haldane proposed to take the degeneracy of entanglement spectrum as the hallmark of topological 
phases, which can be applicable to general topological phases beyond QAFM. 

37 On the other hand, for gapless spin chains, the entanglement entropy diverges as log(L) with L the length of a 
subsystem for which entanglement entropy is defined |57| . 
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B (see Fig. 27). Then, we can express any wave function (state) as [25 

x 

a=l 



(218) 



with non-zero coefficients A Q (> 0). Here, {|a)A} (1 < a < dim^A! Ha'- Hilbert space of A) and 
{|/3)b} (1 < P < dim^B) are orthonormal basis states in the subspaces A and B, respectively: 



(a\a') A = 5 a>a > , (/3|/3') B = Sp^i. 



(219) 



The expansion (218) called the Schmidt decomposition defines the Schmidt coefficients X a and 
describes the entanglement between the two subsystems. The Schmidt number x is the number 
of non-zero Schmidt coefficients and never exceeds the minimum of the dimensions of the Hilbert 
subspaces. From the normalization of \^f), the Schmidt coefficients satisfy X^a=i = 1- The 
"spectrum" of the entanglement energy e a defined by 



(e a > 0) 



(220) 



is called the entanglement spectrum |27| . In terms of the Schmidt coefficients, the (von Neumann) 
entanglement entropy is given by |25| 



Se.e. = - ^2 A " log A « = ^2 €ae ~ 



(221) 



It is interesting to observe that Se.e. m ay be viewed as the ordinary thermodynamic entropy if 
we introduce a fictitious "temperature" T and "partition function" as 



Z(T) = ^V^ e <* , Z(T=1) = 1 



(222) 



and define 



d_ 
df 



{-TlogZ(T)} 



T-s>l 



(223) 







V entanglement cut 








"edge" 




















A 






B 





Figure 27: The VBS state is divided into two parts and the edge degrees of freedom emerge at 
the cut. 
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In the MPS formulation, the derivation of the Schmidt coefficients are rather straightforward. 
Since the MPS is represented as the product of matrices Aj defined on each site, one can easily 
write down an analogue of the Schmidt decomposition (218) (see Sec. 5.2 for the physical meaning 
of the matrix indices): 



D 



(AiA 2 ---A 



N ) aj,, at 



^{AiA 2 ■ ■ ■ Ai) aL , a ■ (Aj+i ■••A. 



N )a,otji ; 



(224) 



a=l 



where D is the size of the matrix Aj and each state-valued matrix A is expanded explicitly in 
terms of the c-numbered matrices A(m) and the orthogonal local basis states |m)j as 



A i = Yl A j( m )\ m ) 



(225) 



m=l 



(d: dimension of local physical Hilbert space). For simplicity of argument, we assume that the 
system is uniform [Aj(m) = A(m)] and defined on a open chain. 

One may think that eq.(224|) already completes the Schmidt decomposition (218) with the 
identification 



\oc)b 



(A1A2 ■ ■ " Ai)ai,, a 

■A N ) 



(A 



(226) 



i+1 



However, the orthonormality condition ( 219| are not always satisfied for the above choice. Nor- 
mally, we use the "gauge ambiguity" [2U [31] to find an appropriate set of edge states in such a 
way that the following overlap matrices may equal to identity: 



£a,(3 — (A^Al ■ ■ ■ A*) ahi a{AlA2 • • ■ A)a L ,/3 
Tta.,13 = (A* + i • • • A* N ) a ,a R {Ai+i ■ ■ ■ ^4jv)/3,a F 



(227) 



This procedure can be carried out both for finite systems and for infinite-size systems by using 
the singular- value decomposition (SVD) [391 SOI Hlj- As a result, we obtain, instead of a single 
matrix A(m), two different MPS matrices AT(m) for the left subsystem A and r(m)A for the 
right subsystem B. The (diagonal) matrix elements of the DxD diagonal matrix A coincide with 
the Schmidt coefficients: 

[A]oa = A Q . (228) 

The MPS characterized by the set of matrices {A, T(m)} (m = 1, . . . ,d) is called canonical 
and automatically completes the Schmidt decomposition (218) as [39],H0lH3] 

D 



I*) = J2 X c c \a) A ®\u)B 



D 



Yl Xa t" " Ar ( m *-i) Ar ( m i)]a L ,a [T(m i+1 )Ar(m i+2 )A ■■■} 

a=l {rrii} 

x |mi)<g>|m2)<g> • • • ®\mi) 



(229) 



a,«R 
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(for a rather comprehensive account of the use of SVD in MPS, see, for instance, |32j). 



In infinite-size systems, the overlap calculation (227) simplifies a lot since C (1Z) reduces essen- 



tially to the left (right) eigenvector V L ^ (VjJ ) with the largest eigenvalue (dominant eigenvector, 



see Fig. 28 for a diagrammatic representation 38 of the transfer matrix Tl (Tr) defined by 

PtiawfiJI = E( Ar *( m ))a^( Ar M)^ ( 230a ) 

rn 

(Tr)«,«;M = E( r *( m ) A )«^ r H A )^ (aJ,a,P = l,...,D) . (230b) 



Specifically, except for the unimportant factors determined solely by the edge states, the overlap 



o — □ □ — □- 



□ — □ □ — □- 



H3 — n □ — □- 



1 LrdJQD — □ □ — □- 



L I 



boundary factor 



Figure 28: (Color online) Norm of MPS on a semi-infinite system. If the MPS is pure, the norm 
is essentially determined only by the dominant eigenvector vj^ (except for the boundary factor 
which depends on the boundary condition imposed on the edge states). Similar relation holds for 
the right semi-infinite system. 



matrices (227) coincide with the (D 2 -dimensional) eigenvectors an d Vd of Tl and Tr, 

respective 



y (see Fig. [28]): 

C &f) oc [Vl 1 ^ , Ke0 oc [V^W . (231) 

If our infinite-size MPS (iMPS) assumes the canonical form, [£]ap oc 5 a p = (l.o)a/3 an d 
[R]ap oc 5 a/3 = (Id) a/3 by definition (with In being the T>-dimensional identity matrix), and 

38 When the MPS is not of the canonical form, we can take T(m) — A(m) and A = 1 (and hence Tl = Tr). 
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therefore the two transfer matrices satisfy 39 [20, 1193] 



D 



a, a a=l m=l 



a 

^2 ^r t (m)A 2 r(m) _ = (lfl)^ (left action) 



(232a) 



and 



m=l 



D d 



E l T *ka;- M (i^W = E E [ r *HA]^ [r(m)A]^ 



0,0 



= 1 m=l 



r t 
|r(m)A 2 rt(m)} 



(232b) 



m=l 



(lo)a,a (right action) 



(see Fig. 29). Notice that these equations can also be regarded as the eigenvalue equations for 



the transfer matrices: Id is the left (right) eigenstate of the transfer matrix Tl (Tr) with the 



eigenvalue 1. Eq.(232) may be thought of as the conditions in order that the iMPS is of the 
canonical form. 



(a) E 



•a a 




■/3 




Figure 29: (Color online) Graphical representation of the two conditions, (232a) and (232b), 
for canonical iMPS. In Ref. |193j . the action of the type (a) ((b)) is denoted by £*(A 2 ) = A 2 
(8(1) = 1). 



On the other hand, when the MPS in question satisfies^] 

C a p OC 8 a fi , TZ a /3 OC 5 a> j3 



(233) 



in the infinite-size limit, the Schmidt decomposition for the infinite chain is obtained very easily 
just by rescaling the MPSs for the subsystems. The normalized state |^) is constructed as 



|*(a,/3)) 



(A1A2 ---A 



N)a,0, 



(234) 



39 By the assumption of pure MPS, the largest eigenvalue is unique |20l 1193] . (MPS is called a pure MPS when 
transfer matrix has non-degenerate maximal eigenvalue. In the infinite-size limit, the pure MPS is reduced to a 
pure state, and hence the name pure MPS.) When the value of the largest eigenvalue is not 1, we can rescale the 
matrices A(m) so that it may be 1. 

40 This is the case for all the (S)MPSs treated in this paper. 
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where A(a,j3) is the magnitude of \^>(a,f3)): 



A(a,p) = \{AiA 2 ---A N ) a ^\ 2 = {A\ ■ ■ ■ A\a\)p >01 {AiA 2 - ■ ■ A N )a,p- 
With these normalization constants, the normalized MPS is written as 

M*L,a R )) = £ J A{aL A f A{a \ aR) m*L,a)) ■ Ma, a R )). 

A(a L ,a R ) 



(235) 



(236) 



Comparing this expansion with (218), one may read off the Schmidt coefficients as 



' A(a L ,a)A(a,a R ) 
A(a L ,a R ) 



(237) 



In the infinite limit, A a is not relevant to the polarization of the edge spins, oil and a R . Therefore 
for the infinite (S)VBS chain, we only need to evaluate the magnitude of the matrix product to 
obtain the Schmidt coefficients, and thus the derivation of Schmidt coefficients is boiled down to 



the computation of the normalization constants, A(a,f3) (235) 



6.3.2 Af = l 

To substantiate the topological stability of the SVBS chain, we investigate two type-I chains with 
distinct bulk superspins, 5 = 1 and 5 = 2. 
(i) Superspin 5 = 1 

From the formula (237), the entanglement spectrum of the 5 = 1 type-I state (on an infinite 
chain) is derived as 



Ab = Ai 



\ 2 
Af 



A, 2 



A 2 



1 



1 3 

4 4V9 + 8r 2 ' 
3 



2 2V9 + 8r 2 ' 



(238a) 
(238b) 



They are shown in Fig. 30 (left-figure) with the entanglement entropy (right-figure). At r = 0, 
the SVBS state reduces to the 5 = 1 VBS state and reproduces both entanglement spectra and 
entanglement entropy of the 5 = 1 VBS chain, i.e., Ab 2 — > 1/2 and See — > ffi2, which should be 
compared with (217) for 5=1. Similarly, in the limit r — > oo, the SVBS chain reduces to the 



MG chains, and the entanglement entropy of SVBS chain also reproduces the finite entanglement 



entropy of the MG chain (the right-figure of Fig. 30). 



As in the entanglement spectra of the left of Fig. 30 we have two distinct entanglement spectra, 
one of which is doubly degenerate spectrum (blue curve) for the "bosonic" Schmidt coefficients 
corresponding to those of the original 5=1 VBS chain (238a) and the other is the non-degenerate 



spectrum (red curve) for the "fermionic" Schmidt coefficient (238b). The existence of such two 



types of entanglement spectra is a salient feature of the SUSY state and can be readily understood 
based on the following edge state picture. For 5 = 1 SVBS chain, its edge superspin states are 
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Figure 30: The change of the entanglement spectra as a function of r (left panel) and entanglement 
entropy (right panel) of the 5 = 1 type-I chain. The insets are the results for the bosonic-pair 
VBS chain. (Figure and Caption are taken from Ref,|10|.) 



given by the UOSp(l\2) multiplet with superspin S dge 
states with spin, 1/2 

(Sedge = 1/2) ffi (Sedge = 0) 



1/2 that consists of the ordinary SU (2) 



(Fig. |3Tp : 

(Sedge = 1/2) 



or 



(239) 



(240) 



The 5 = 1/2 SU(2) edge spin states generate the double degeneracy in the entanglement spectra, 
while the 5 = SU(2) edge spin state gives the non-degenerate one. Due to the existence of 1/2 
spin edge degrees of freedom, the double degeneracy is guaranteed (we will give a detail discussion 



in Sec 6.4), and hence we find that the type-I SVBS chain is in the topological phase. 




Edge superspin 
BJ S ed 9 e = 1/2 




Bulk superspin 

S = l 



Figure 31: (Color online) There always exit half-integer and interger edge spin states as the 
super-partner of the SUSY. The half-integer edge spin states play a crucial role in the stability of 
topological phases. 



To highlight the effect of SUSY, we consider the non-SUSY (i.e. purely bosonic) cases and 
replace the fermion operator /, with the boson operator q ([cj,cj-] = <%) to define the following 
boson-pair VBS chain, 

l*b. P .> = U^l+i - 6 H+i - rc M+i)|vac>. (241) 
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0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



1 + r 2 



Figure 32: The behaviors of the Schmidt coefficients (243) and the entanglement entropy (inset) 
of the 5 = 2 type-I SVBS chain. (Figure and Caption are taken from Ref. [10].) 



As the insets of Fig. 30, we depicted the entanglement spectra and entanglement entropy. The 
crucial difference to the SUSY case will be apparent in the limit r — > oo, l^b.p., 
simple product state, 



is reduced to a 



l*b.t 



n 



cjlvac) 



(242) 



and the entanglement entropy vanishes (the inset of the right figure of Fig. 30). 
(ii) Superspin 5 = 2 

Next, we examine the entanglement spectrum of the 5 = 2 type-I SVBS chain. Though the 
topological phase of its bosonic counterpart, 5 = 2 VBS chain, is fragile under perturbation |190j . 
the 5 = 2 type-I SVBS chain itself is topologically stable with finite amount of hole doping. As 
we shall see below, SUSY plays a crucial role for the stability of the topological phase. At r = 0, 
5 = 2 SVBS chain is reduced to the 5 = 2 VBS chain, while in the limit r — > oo the SVBS chain 
is reduced to the partially dimerized chain. We have the following five Schmidt coefficients for the 
5 = 2 SVBS chain: 



A2 2 — A3 2 — - + 

6 



Ab 2 = A1 2 

Af 2 = A4 2 = A5 2 = - — 



1 5(4 + V25 + 24r 2 ) 

6(25 + 24r 2 + 4^/25 + 24r 2 ) ' 
5(4 + ^25 + 24r 2 ) 



4 4(25 + 24r 2 + 4^25 + 24r 2 ) ' 



(243a) 
(243b) 



Thus, the five Schmidt coefficients are split into the triply degenerate (243a) and doubly degen- 



erate (243b) spectra showing distinct behaviors in Fig. 32 Again, such splitting of the Schmidt 
coefficients are readily understood by the edge state picture for the SUSY chain. For the 5 = 2 
type-I SVBS chain, the edge superspin is given by 5 e d ge = 1 that consists of the SU (2) edge spin 
5 e dge = 1 and 5 edg e = 1/2 [Fig. [33]: 



edge 



1) 



edge 



1) © (5 edge = 1/2) 



(244) 
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or 



5 = 30 2. 



(245) 



The /Sedge = 1 degrees of freedom generate the triple degeneracy in the entanglement spectra, 
while the Sedge = 1/2 degrees of freedom give the double degeneracy. Due to the existence of 
SUSY, the 5 = 2 SVBS chain necessarily contains the S e d ge = 1/2 edge spin degrees of freedom 
that do not originally exist in the S = 2 VBS chain [Fig. 33 , and they guarantee the double 



degeneracy in the entanglement spectra, i.e. the stability of the topological phase [61] [62] . 




Edge superspin 



s, 



edge 



i 



Bulk superspin 
5 = 2 




Figure 33: (Color online) S e d ge = 1/2 generates the stability of topological phase of the 5 = 2 
type II SVBS chain. 



From the above demonstrations for 5 = 1 and 5 = 2 type I SVBS chains, one may see that 
regardless of the parity of bulk superspin, the SUSY introduces the half-integer edge-spin states 
that necessitate at least double degeneracy in the entanglement spectrum. 



6.3.3 TV = 2 

It is also straightforward to calculate the Schmidt coefficients for the type-II SVBS chain: 



Ai 2 — A2 2 



1 



4 4V9 + 10r 2 + r 4 ' 



A, 



A, 



1 



+ 



4 4V9 + 10r 2 + r 4 ' 



(246) 



The four Schmidt coefficients are split into two groups showing distinct behaviors [Fig. 34 
cording to the SU{2) decomposition of the UOSp(l\2) edge superspin state: 



ac- 



(5 ed ge = 1/2) = (Sedge = 1/2) © (Sedge = °) © ( S edge = 0), 



(247) 



or 



(248) 



The first 2 on the right-hand side of (248) corresponds to the doubly degenerate blue curve in 



Fig. 34, while the remaining 1 © 1(= 2) represents the doubly degenerate red curve. 



66 



r = v / 3 : 

I 

0.0 0.2 0.4 0.6 0.8 1.0 



1 + r 2 

Figure 34: The behaviors of the Schmidt coefficients and the entanglement entropy (inset) of the 
S = 1 type-II SVBS chain. 



6.4 Supersymmetry-Protected Topological Order 

In this section, we clarify the relation between the structure of the entanglement spectra discussed 
in the previous sections and the symmetry- protected topological order [60J. To this end, we use the 
MPS representation established in Sec. [5j Because of its simplicity in the entanglement structure 
and the wide applicability to gapped states in ID, the MPS approach provides us with a powerful 
tool in investigating the topological phases in ID [611 [621 [6H |65l 1551 155]. 



6.4.1 Symmetry Operation and MPS 

For later convenience, we give a quick summary of some tools used in the MPS approach. As has 
been discussed in Sec. [5] and 6.3.1, any MPS on an open chain may be written as 



|MPS) = (g) Ai = • • • <g> A-i <g> Ao <8> Ax <S> A 2 <8> • • • , (249) 

i=— oo 

or in terms of the c-number matrices A(m) [Aj = Ylm=i A?( m )l m ).?> ^ being the dimension of the 
local physical Hilbert space], as 

|MPS) = ^2 {" " " A{m- 1 )A{m )A{m 1 )A{m 2 ) • • • } |mi)<g>|m 2 )<g> • • • ®\m L ) (250) 

{mi} 



As has been mentioned in Sec. 6.3.1, we can "gauge-transform" this into the canonical form 
which conforms with the Schmidt decomposition; on the left of the entanglement cut, we use 
the MPS matrix A-^(m) = AT(m), and on the right A-^{m) = r(m)A. In the following sections, 
our arguments will be based on the canonical form characterized by the MPS data {A, T(m)} 
(to = 1, . . . , d). 
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Now let us consider unitary operations on MPS. A given MPS (|^mps)) is sai d to be invariant 
under a unitary operation if the product of local unitary operators <g>u leaves the MPS invariant 
(up to an overall phase) |193| 1194] : 



u ® u <g> • • • ® u|#mps) = e lQg - s - |*mps) 
The local unitary operation u acts on A in a site-wise manner: 



m)\m)i 



E A ) ^2(m\u\n)iT(n) > |m)j , 

m \ n ) 



(251) 



(252) 



where we have used the completeness relation on each site, ^ m \m)i{m\i = 1. If the symmetry 
operation is anti-unitary (like time-reversal), the complex-conjugation T(m) h-> T*(m) should be 



taken in eq.(252) 



Then, it can be shown that the above unitary invariance is equivalent to the existence of the 
following unitary operator U acting on T(m) [193J: 



U Q$l U i 



«|*MFS> = e^-I^MPs} 



^2(m\u\n)T(n) = e Wu U ] T{m)U , (253) 



n=l 



where 8 U denotes a n-dependent phase. Once the unitary u is given, U is uniquely determined 
[193J up to an overall phase. To be more precise, if we define a generalized transfer matrix 



[T (u) Ua;/M= E [Ar*(m)]^[Ar(n)] ai/3 ( 



m\u\n) 



(254) 



m,n=l 



the Z)-dimensional unitary matrix U is given essentially by its left eigenvector with the largest 



eigenvalue 41 e ldu (see Eq.(270) and the explanations around it). Since U leaves the MPS |^mps) 
invariant, it is natural to assume that U does not change the physical entanglement spectrum A, 
i.e. 

[U , A] = . (255) 



Physically, the above relation ( 253| ) implies that the original symmetry operation (acting on the 
physical Hilbert space at each site) "fractionalizes" into the ones (U and IP) which act on the 



edge states on both ends of the system. The equation ( 253 ) plays a crucial role in the following 
discussions. 



6.4.2 Case of SMPS 

Now we extend the arguments developed by Pollmann et al. |62l I6T] for the bosonic MPS to the 
SUSY case. First remember that degeneracy in energy spectra of quantum mechanical Hamilto- 
nian can be attributed to some symmetry of the Hamiltonian. Meanwhile, since the entanglement 
spectrum is solely determined by a ground-state wave function, the degeneracy of entanglement 

41 We have assumed that T(m) is rescaled in such a way that the largest (or, dominant) right-eigenvalue of Tl,r 
is unity. 
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spectrum is expected to stem from some symmetry of the ground-state wave function. Indeed, 
several discrete symmetries are identified to guarantee degeneracy in the entanglement spectrum. 
Since the degeneracy of the entanglement spectrum is the hallmark of the topological order, it 
is said that the topological order is protected by the symmetry of the ground-state wave func- 
tion, and hence the name symmetry protected topological order. In the following, we will show 
that the SUSY guarantees the existence of at least two-fold degeneracy of entanglement spectrum 
regardless of the parity of superspin, provided that at least one of the three symmetries, inversion, 
time-reversal and Z 2 X Z2, is present. We will also find that the (S)MPS formalism plays a crucial 
role in the discussions. 



As has been discussed in Sec 5.2, the SMPS is generally represented as [9] 

\V)=AiAi'--A L , (256) 

where A% {i = 1, 2, • • • , L) are supermatrices defined on the sites i. Then, we follow the same steps 
as in Sec. |6.4.1 to transform the SMPS into the canonical form. When a given unitary operation 



leaves the MPS invariant [in the sense of eq.(251)], the unitary operation u fractionalizes into U 
and U\ and acts like 

d 

J2(m\u\n)T(n) = e idu U j T(m)U . (257) 

n=l 

In fact, the local symmetry operation need not be unitary (as will be in the following sections). 
In this case, the left-hand side may be replaced with the general form T'(m). Therefore, the most 
general form reads as 

r'(m) = e i8u U ] T(m)U , (258) 

where the phase Ojj depends on the symmetry operation considered. In the above equations, m 
labels both bosonic and fermionic states, i.e. m = i,a, and T(m) are given, for UOSp(l\2), by 





= x,y,z) 

(259) 

(a = 0i,0 2 ) , 



Mi(i) 
M 2 (i 

Ni(a 
N 2 (a) 

where M\,M 2 ,Ni and N 2 are c-number matrices. 

Here, a remark is in order about the form of U. The (c-number) unitary matrix U in (257) 
may be postulated as: 

where U\ and U 2 are unitary matrices that act on the two bosonic subspaces having different 
fermion numbers. The reason for choosing the above form may be seen as follows. First we note 



that eq.(257) implies that the MPS on a periodic chain transforms like 

^ STr(uU 1 A 2 ■ ■■Azn+xU) , (261) 
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where the supertrace STr is defined for general supermatrices as [see also eq.(176b)] 

,(2) 



STr 




TrA 



(i) 



TrA\ 



(262) 



We expect that for a periodic chain, which does not have edges, is fully invariant under the unitary 
operation, i.e., the expression (261) coincides with the original MPS (up to an overall phase). 
While in the case of bosonic MPS, this, combined with Tr(AB) = Ti(BA), immediately implies 
®u\*I>) cx |vl/), the relation STr (AB) = STr(BA) holds only when A and B are super-matrices 
(that contain the Grassmann-odd blocks in their off-diagonal parts): 



|*) ^4 STr(rtAiA 2 ■ ■■A 2n +iU) = STr(^!^ 2 • ■ ■ A 2n+ iUU^ . 



(263) 



In fact, if an arbitrary pair of two super matrices A and B were merely c-number matrices, A 
and B, in general, would not commute inside STr(-): STr(AB) ^ STr(BA). In order to satisfy 
STr(AB) = STt(BA) only with c-number matrices, either A or B is forbidden to have c-number 
components in the off-diagonal blocks. 

For later convenience, we derive a useful property of pure canonical MPSs [24, 44] . In the 
following sections, we assume that the MPS in question is defined on an infinite-size system. The 
equation (258) is the most general statement about how a given symmetry of MPS is realized 
by a projective representation. However, sometimes it happens that T'(m) = T(m) for a certain 
symmetry operation and in these cases we can draw an interesting conclusion about the properties 
of U [5T 1 IS2" | [TU]. Suppose that we have a pure MPS whose canonical form is characterized by the 
MPS data |24} 144] {A, T} and that it satisfies the following relation for some unitary matrix U: 



T(m) 



je. 



u U ] T{m)U . 



Since the MPS is canonical, the following holds [see eq. ( [232a] ]: 

^r t (m)A 2 r(m) = 1 D . 



Physically, it states that the Z) 2 -dimensional vector Vt 



a, b 



>ab 



(l<a,b<D) 



is the dominant left-eigenvector of the left transfer matrix (see Fig. 29 ) 

m 

Plugging rt(m) = e- i6u U^(m)U into ( |265| , we obtain: 

-&v 



or equivalently 



YU ] T ] {m)U h 2 T(m) = 1 D 

m 

J2r\m)AUAT(r 



TO 



u 



(264) 



(265) 



(266) 



(267) 



(268) 



(269) 
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This implies that the DxD unitary matrix 

(U)u = £ {l®U} aa . u = Y, KlVU , (270a) 

a a 

when viewed as a D 2 -dimensional vector, is the left-eigenvector of Tl with the eigenvalue e i(>u : 

UT L = e i6u U. (270b) 

Since, by assumption of canonical MPS, Id is the only left-eigenvector having the eigenvalue 
|A| = 1, we conclude 

e i8 u = i > u = j* a l D . (271) 
Since, in deriving the above, we have only assumed that the (infinite-system) MPS in question 



is pure and takes the canonical form, ( 271 ) holds for any MPS (including SMPS) satisfying the 
assumption. 

6.4.3 Inversion Symmetry 

Now let us look at what inversion symmetry I with respect to a given link implies about the 
structure of the entanglement spectrum. It is convenient to consider an inversion transformation 
on a circle: 

y = STr(AiA 2 ---A 2 n + l). (272) 
The inversion on a given link transforms the SMPS chain as 

m = STr(^ 2n+1 ^ 2 „ • • -Ai). (273) 



By the relation STr(MiM 2 ) = STr((MiM 2 ) st ) = STr(M| t M 1 st ), ( |273| ) can be rewritten as 

m = STv(AfA s 2 t ---A s 2 t n+1 ), (274) 
where 'st' stands for supertranposition defined by 

(275) 



Mi Ni\ _ ( Ml N% 



N 2 M 2 \-Nl M\ 



It is easy to verify that applying the supertrace twice does not return a supermatrix to the original 
one: 

\\N 2 M 2 I I -N 2 M 2 \N 2 M 2 y 1 



where the matrix 



P^( U ° ) (277) 



-1 



2 , 



has been defined in such a way that its adjoint action P{-)P multiples the fermionic blocks by a 
factor (—1). [li and 1 2 correspond to the unit matrices of two bosonic subspaces and should not 
be confused with the DxD unit matrix l^-l 
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Therefore, at the level of the local MPS matrix, the inversion X is realized as 



(278) 



Finally, from eq.(258), one sees that when \& has the inversion symmetry, the .A-matrix should 



satisfy the relation [10 



r(m) st = e idl U\T{m)U I . 



(279) 



Note that at this stage, we can not apply the argument in Sec. |6.4.2 since eq.(279) does not assume 
exactly the same form as eq.(264) [the LHS is not Y{m)]. 



In order to obtain the relation of the form (264), we combine eq.(279) with the fact that the 
link-inversion squared to unity X 2 = 1. Applying supertransposition (-) st to (279) once again and 
using (A st ) st = PAP [eq.d276b], we obtairp 



r(m) = e 2idl (C/ / PC/;) t r(m) (UiPUj] 



This is of the form of eq.(264) and then (271) immediately implies that 



(U I PUf)=e t *'l D , e 



l^e 4 ' 



±1 



After multiplying Uj from the right and making transposition, we deduce 



Uj 



PUj y 



P l Ui 



Ui 44> e" 



±1 



[note PUi = UiP by eqs.fl260[) and ( |277| ]. Therefore, we obtain {10 

U\ = ±PU! . 



(280) 



(281) 



(282) 



(283) 



Eq. (283) states that when U\ is symmetric (anti-symmetric), U2 is anti-symmetric (symmetric). 



It should be noted that unlike in the bosonic case |6H 162] , the symmetry constraint is imposed on 
each of the spin-S" "bosonic" sector U\ and 5 — 1/2 "fermionic" sector U2 in the case of SUSY. 

When U\ is anti-symmetric, for instance, the sector with (—1) F = +1 must have a special 
structure in its entanglement spectrum. In fact, by computing the determinant of U\ 



det U x = det U{ = det(-[/i) = (-l) dl det Ui (di : dimension of Ui 



(284) 



one can immediately see that d\ should be even (the same argument applies to U2 as well). From 
this, one can conclude that either fermionic (when the sign + occurs) or bosonic (— ) sector has 
even-fold degeneracy in each entanglement level, which we can use as the fingerprint [601 W\\ [62] 
of the SUSY-protected topological order [TO] . 



2 In the bosonic case, the matrix P is not necessary |61l 162] . 
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6.4.4 Time-Reversal Symmetry 



We can draw a similar conclusion for the time-reversal symmetry T ■ Let us recall the time-reversal 
operation for the superspin, Si and S a 



St ^ (e^KjStiKe-^y) 
S a ^ {e inS yK)S a (Ke- i7rS v 



Si, 



(285) 



with complex conjugation operator K. The set of equations (285) implies that, for integer super- 
spins S, T = Ke~ mSy satisfies^ 

r 2 = r = (-i) F . (286) 

where the d-dimensional matrix V, which acts on the physical Hilbert space and multiples a minus 
sign when the state is fermionic (i.e. when the fermion number F = odd), is analogous to the 
D- dimensional matrix P [see eq.(277|] acting on the auxiliary space. 

From (285), one sees that, in terms of the canonical matrix, the time reversal transformation 

r 



is represented as 



r( m ) A r(m)' = Y,{ Ry (^UnT(ny, (287) 

n 

where m,n = i,a and R y (ir) = e mSy . For instance, R y (ir) is given, for superspin S = 1, by 



R y (n) 



(-1 











o\ 





1 

















-1 




















1 


\o 








-1 


0/ 



(288) 



Again by eq.(258), the time- reversal invariance implies that there exists a unitary matrix Ut that 
satisfies 

T. 



r(m) A ^i^»r*(n) = e*4r(m)f/ T 

n 

Applying this twice and using T 2 = V, we obtain 



(289) 



(V)uT(l) = J2l RV ^)]lm J2[R y (-K)}mnT* 



, n 



m=l 

d 



.n=l 



m=l 



(290) 



e-*[4 0][E»W] [m r(m) ^ 



, m=l 



= {U t Ut} ] T(1){U t Ut} 

r(0 = {UtUtV (P)a r(0 {tW} • 



43 Recall that, in the ST/ (2) case, T 2 = +1 for integer S. When superspin S is half-odd-integer, the time reversal- 
operator satisfies T 2 = —P. (This is indeed a generalization of T 2 = — 1 for the SU(2) half-integer spin case.) 



73 



Using the relation (V)uT(l) = PT(l)P, we may rewrite (290) into the following form 

T(l) = {U t PUtV r(J) {U t PUt} , (291) 



to which we can apply the argument presented in Sec. |6.4.2 [see eq.(264)]. Thus we obtain 

U T PU£ = e i$T l . (292) 



(the value of 0t can not be determined by this method). This is of exactly the same form as (281 ) 
and we deduce the same conclusion: 

C4 = ±PU T , e iq>T = ±1 • (293) 



Therefore, as in the previous case [see eq.(283)], we see that time-reversal symmetry guarantees 
the existence of at least double degeneracy in the entanglement spectrum of (either of the bosonic 
or the fermionic sector). 

6.4.5 Z 2 x Z 2 Symmetry 

Lastly, we consider the two independent tt rotations around the x [u x (tt)] an d the z axes [u z (ir^ ' 1 
As will be seen below, u x {n) or u z (tt) alone does not lead to any significant conclusion. However, 
the combination of the two |6HI62j leads to a similar conclusion about the entanglement spectrum. 
Under the tt rotation around the x (z) axis u x (tt) (u z (tt)), SMPS transforms as 

r( m ) r(m)' = ^ R a mn ^)T{n) (a = x,z). (294) 

n 

When the SMPS respects such a symmetry, we have 

Y,[R a ^)]mnT(n) = e^UlT(m)U a , (295) 

n 

where (R a (ir)) 2 = V for integer superspin S. Again we follow the same steps as before to show 

e 2i8 a = 1 ^ e i9 a = ±1 

(296) 

U a PU a = e t9 'l (a = x,z). 

Note that the phase factor e** a can be absorbed in the definition of U a and hereafter we assume 
Ua = PU a (a = x,z). Unlike in the previous cases, this relation alone does not give any useful 
information about U a . 

On the other hand, the combination of the "commutation relation" u x (tt)u z (it) = 7 ? M 2 (7r)« a; (7r) 



and eq.(295) implies, in terms of r(m), 

r(m) = (U x U z UtPUi)T(m)(U z PU x Uiul) (297) 

and hence gives 

{U z PU x )(UtuD =e**«l . (298) 



J Note that u v (tv) is redundant since u x (tt)u z (tt) — u z (tt)u x (iy) — u u (ir). 
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In obtaining eq.(297), the phase factor e lSxz = e i x e i " z appears. However, it cancels out in the 



final expression (297). 

Since the phases of U x and U z have been already fixed, the phase of U X U Z cannot be arbitrary 
and has a definite physical meaning. By multiplying U Z P and U x from the left and the right of 



(298), respectively, and using Ul = PU a = U a P repeatedly, one can show 

U x U z = e^U z PU x , (299) 



which is combined with (298) to give e l ® xz = dbl. To summarize, the two unitary matrices U x and 
U z obey the following commutation relation: 

U X U Z = ±PU Z U X . (300) 



In terms of the block components U a \ and U a2 , eq.(300) reads as 



U Xtl U Zll = ±U zA U xA , U X , 2 U Z>2 = TU Z , 2 U X>2 . (301) 



As in the previous cases, (301) immediately implies at least double degeneracy of each entangle- 



ment levels in the sector (1 or 2) taking the minus sign in the above equation. 
6.4.6 String Order Parameters and Entanglement Spectrum 

In the previous sections, we have shown that, as in the purely bosonic cases |60| [6T| [62], such 
elementary symmetries as inversion and time-reversal protect the Haldane phase from collapsing 
into trivial gapped phases on the basis of the assumption that the even-fold degeneracy in the 
entanglement spectrum either in the fermionic or in the bosonic sector is the entanglement finger- 
print of the topological "Haldane phase" in SUSY systems. On the other hand, the string order 



parameters [see Sec. 6.1 have been used traditionally to characterize the Haldane phase [IS! H9J. 
Then, a natural question arises: is there any connection between the description by the string 
order parameters and the modern characterization in terms of the entanglement spectrum? In 
this section, we give the answer to this question. 

Let us first consider the structure of the string order parameters and Coring [eq.( |206 )] 



from the MPS viewpoint [23, 193J. When we evaluate them using MPSs, we encounter the following 
matrices |23j : 



[T a h >a .- M ^ £ [A*(m)]^[A(n)} a ^(m\S a \n) (a = x,z) 

m,n=l 

d 

[Tstr ms ]aM= £ [A*(m)]^[A(n)] a ^(m\exp(i7TS a )\n) (302) 

m,n=l 
d 

[T^in g ]aM= E [A*(m)]^[A(n)] a ^(m\S a exp(i7rS a )\n) (a = x,z) 

m,n=l 

as well as the usual transfer matrix T. By using these matrices, the MPS expression of the string 
order parameter Oaring ^ s §i ven as 

T^T^T^^T* T Nr , (303) 
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where we have omitted the denominator necessary to normalize the MPS. The parts T Nh and 
T Nr are easy; for the canonical MPS, they reduces, in the infinite-size limit, to (see Sec. 6.3.1): 



[T Nh ] 



ol,ol;/9,/3 



IT** WflaA • (304) 



The boundary dependent factors 5q LjQl and #a R)l g R are cancelled by those coming from the de- 
nominator. Therefore, the quantity which we have to compute is 



T z (Ta. • ^-^T 2 

A string \ A string; J- 



a,a;P,P 



(305) 



Since we are interested in the long-distance limit \i — j\ /* oo, we need know the asymptotic 
behavior of the string (String) ■ First we note that T str i ng may be thought of as the overlap 



(MPS|u<g>u®---®u|* 



MPS/ 



(306) 



with u = exp(iirS a ) (a = x,z) [see eq.(251)]. Then, it can be shown |193j that in order for the 
string (T str i ng )' l-J ' not to vanish in the long-distance limit, the MPS should be invariant under 
both of the n- rotations u x (^) and u z {tt). 

Then, from the discussion in Sec |6.4.5 (we just replace P i— >■ 1 to obtain the results for the 



purely bosonic case) , we see that there exists a pair of unitary matrices satisfying ( 295 ) and 



UM Z = ±U Z U X 



(307) 



see eq.(300)]. Unfortunately, we cannot tell whether the even- fold degeneracy in the entanglement 



spectrum, which is the fingerprint of the topological Haldane phase, happens or not since we do 
not know which of the ± signs is chosen. 

Now we show that when the string order parameters are non- vanishing Coring ^ ^> the 



minus sign in fact realizes (i.e. U x and U z anti-commute) in eq.(307) and the entanglement 



spectrum exhibits the degenerate structure. To this end, we investigate eq.(305). Since the 
largest (right) eigenvalue of T str i n g is e l9a [see Sec. 6.4.1 , (Turing)! 1- -?! reduces essentially to a phase 



( e ifl a )|i-j| = (±i)[*-il. The price to pay is the following boundary factors appearing at the two 



end points of the string correlators (see Fig. 35 ) 



D 2 



{ ^string I E ^R,n^L,n | (^strim- . 



\i-j\rpz 



a,P 



,n=l 



a,a;P,p 



\i-j\/-oo 



E {(String vg)(vgr 2 )} = E {(T* tI - m Ji®ut\ 1)(1 {l®U z } T Z )\ 

1 — ' l- ) a,a;P,p * — ' V V ) ) a,a; 

a,P a,P 



;P,P 



(308) 



where V^] (VjJ*^) denotes the left (right) dominant eigenvector of T str i n g and the tensor notations 
are defined as 



[l® B ha;P,P ~- 



5 a,p \ B \aP ' 



[{1®B} l] &ja = E 5^ [B] ap b u = [B% a . ( 309 ) 

P,P 
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In fact, the eigenvalues of T s t r i ng precisely coincides with those of T for all the spin-S 1 VBS states 
and the different behaviors in the string order for the odd-S 1 and the even-S* chains comes only 
from these boundary factors. 

To see whether the boundary factors are non- vanishing or not, we consider the right-boundary 
factor (1 {l(g>U z } T z ) of 0g tring . First we rewrite it by using (see the second figure of Fig 36): 



S z 



(u\.u x )S z (u\.u x ) = u x {-S z )u x (u x 



(310) 



The unitary operators u x and u x appearing on both sides of — S z can be absorbed into the MPS 



matrices by using eq.(295) (the third figure of Fig 36). By re-arranging the unitary matrices U X U : 



with the help of eq.(307) (the fourth figure of Fig. 36), we arrive at the expression: 



1{1®U Z }T Z = 1 [l®{U x U z Ut)} (~T Z ) 

= 1 {l®(±tf z l7 x tft)J (-T z ) = . 



(311) 



Therefore, we see that the boundary factor (1 {1®U Z } T z ) vanishes (and so does the string order 
parameter 0^!°° ) when U x and U z are commuting (i.e. when the minus sign in the last expression 
is chosen). By the explicit construction of U x and U z , we can easily see that for the even-S VBS 
state, U X U Z = +U Z U X holds. Therefore, it immediately results that the string order parameters 
identically vanish for the even-S VBS state solely for a symmetry reason. 



S a c" 



9f 
1 1 1 



rpa 



19 

: ' ^ 

•fa 
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ill: 




si 





td 



Figure 35: (Color online) Diagrammatic representation of the main part of string correlation 
function j (^ s tringVRi)(VL U j T z )j. vj^j denotes the dominant eigenvector of T str i ng . (Figure 
has been adapted from Ref. |10j.) 

In short, the existence of non- vanishing string order parameters Coring ^ an d C'string ^ ^ 
implies that the existence of the two anti- commuting unitary matrices U x and U z , and thereby 
guarantees the even-fold degeneracy in the entanglement spectrum. To put it another way, the 
string order parameters work as the sufficient condition for the topological Haldane phase. It is 
crucial that both of the string order parameters are finite. For instance, if we "deform" the original 
St/(2)-invariant VBS state by using the quantum group U q (su(2)), we obtain a VBS state with 
uniaxial anisotropy [22J, where one of the string order parameters (0^^° ^ 0) vanishes while 
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±u z u x 



Figure 36: (Color online) Rewriting the boundary factor 1{1®U Z }T Z (for a = z) using u x [see 
eq.(311)]. When U x and U z are anti-commuting, the minus sign coming from u x S z u x = —S z is 



cancelled and an overall plus sign is recovered. (Figure has been adapted from Ref. [TO].) 



the other is still finite [116] . In this case, the unitary U x does not exist and the existence of 
the degenerate structure in the entanglement spectrum is no longer guaranteed. In fact, explicit 
calculation shows that the two-fold degenerate entanglement levels in the S = 1 VBS state split 
into two non-degenerate levels in the deformed state. 

It is straightforward to generalize [10J this to SUSY cases by taking into account the appearance 
of the P-matrix. Again by using u x (7t)S z u x (tt) = —S z and U X U Z = ±PU Z U X , the second factor 



of the right-hand side of eq.(308) can be recasted as: 

1{1®U Z }T Z = i{i®(cy/*tft)}(-r*) 

= l{l®(±PU z U x Ul)} (-T z ) = ^1{1®PU Z }T Z . 



(312) 



Since this implies 

£ {1 {1®^} T Z } M = Y, {{1*U M } T*} aM 

P a,P 



■o- 



■u- 



aeB 



-o- 



u. 



Vz 



-o- 



(313) 



-o- 



-L> 



■O- 



we see that one of the two components (bosonic or fermionic) vanishes just by symmetry: 

( > □ v 





f 








Q 




S'\ 


i 








r 


h 


) 



S* 



=8= 



when e % ® xz = +1 



(314) 



-D- 



when e l ® xz 



-1 . 
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In conclusion, we have established the connection between the string order parameters, which 
have been commonly used [195\ 1162] 11961 1197] [T98J to characterize the Haldane phase, and the 
entanglement spectrum [27J which is a modern tool to look at the topological properties in the 
bulk. As has been shown in Sees. 6.4.3 6.4.5 SUSY guarantees, regardless of the parity of the 
bulk superspin S, the existence of the even- fold degeneracy in the entanglement spectrum of either 
of the bosonic or the fermionic sector, provided at least one of the three symmetries (inversion, 
time-reversal and Z2 XZ2) is present. Therefore, in contrast to the usual bosonic (i.e. SU{2)) VBS 
states, in which the Haldane phase is stable only for odd-integer spins, in the SVBS states, SUSY 
plays the crucial role in protecting the topological Haldane phase. The "revival" of the string 
order parameters upon doping the system (discussed in Sec. 6.2) may be naturally understood by 
the above argument. 

As has been seen above, that both of the two string order parameters are non-vanishing is 
the sufficient condition for the symmetry-protected Haldane phase. However, one may seek more 
faithful order parameters and, in fact, some results have been obtained in this direction. See, for 
instance, Refs. [199J and [200] for recent attempts at finding better order parameters. 



7 Higher Symmetric Generalizations 

In this section, we extend the previous SUSY formulation to the higher symmetric UOSp(l\A) 
case. We begin with the construction of fuzzy four-supersphere based on the UOSp(l\4) algebra 
|123j . As the coordinates of fuzzy four-sphere correspond to the SO (5) gamma matrices [92J, the 
coordinates of the fuzzy four-supersphere are constructed from "super gamma matrices" of the 
UOSp(l\A) algebra. Next, based on the UOSp(l\A) structure, we derive the UOSp(l\A) SVBS 
state [TO] whose bosonic counterpart is the SO (5) VBS state [Ml EES! [TO90 We also develop 
the SMPS formalism for the UOSp(l\4) SVBS state and investigate the topological properties. 

7.1 Fuzzy Four-supersphere 

Here, we utilize UOSp(l\4) algebra to construct fuzzy four-supersphere with M = 1 SUSY. It may 
be worthwhile to first point out a nice correspondence between algebras and fuzzy spheres. For 
fuzzy two-sphere, we have 

SU{2) ~ USp(2) -»■ 5 f 2 , UOSp{\\2) -> sf 2 , (315) 
and for fuzzy four-sphere, 

50(5) ~ USp(A) 5 f 4 , UOSp(l\4) -> sf {2 . (316) 

The UOSp(l\4) algebra is constituted of fourteen generators: 

dim[«o«p(l|4)] = 10|4 = 14, (317) 

45 There still exist the correspondences among fuzzy geometry, QHE, and VBS in the higher dimension case. The 
SO(5) VBS corresponds to the Laughlin-Haldane wave function in 4D QHE [75] 1201] . and 4D QHE realizes the 
fuzzy geometry of fuzzy four-sphere. 
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ten of which are the bosonic generators T ab = —T ba (a,b = 1, 2, ■ ■ ■ ,5) (the SO (5) generators) and 
the remaining four are the fermionic ones T a (a = 1,2,3,4) (the SO (5) spinor). They amount to 
satisfy the following algebra: 

[r a &, T c d] = i(5 ac T bd - 5 ad T bc - 5 bc T ad + S bd T ac ), 

{r Q ,r^} = ^(c 7afe ) a/3 r ab , (318) 

a<b 

where C is the SO(5)(~ USp{4)) charge conjugation matrix 

to 

and j ab are the 50(5) matrices. In the following discussion, we take the 50(5) matrices as 

1 Ar 3 \ 1 ( (72 \ 1 (a ± \ 

712 =2^0 a 3 J' 713 =2^0 -a 2 J' 714 = 2 ^ - ffl J ' 

1 / -(tA l/ffi o\ 1 /<r 2 \ 

715 =2^ J' 723 = 2 ^0 aj' 724 = 2 ^0 -<x 2 J ' 

1 / -<r 2 \ 1 (a 3 \ 1 / -<7 3 \ 

725 = 2 ^-a 2 J' 734 = 2 ^0 -osj' 735 = 2 ^-a 3 J' 

The f705p(l|4) quadratic Casimir is given by 

£ = y~] r a6 r ab + c a( g r Q r0. (321) 

The fundamental 5-dimensional representation matrices of uosp(l|4) are constructed as follows. 
From the (bosonic) "gamma matrices" of UOSp(l\4) algebra 

r *=( 7 n a °X ( 322 ) 



with 
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11= [-ia 1 ' 72 = --2 
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we can derive the SO(5) generators 



r ab = -i-[T a ,v b ]= r a Q b ° 



(324) 



V2\-(cr a y or 



The fermionic "gamma matrices" are also constructed as 

r, 

where 

n 



They satisfy the "hermiticity" condition 

r| = r a , r* b = r ab , r*, = c a pTp, 



(325) 



(1\ 




AA 




M 
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, 73 = 





, n = 
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w 




W 



(326) 



(327) 



4|2 

the coordinates of S { , X a and Q , as 

X = —^T ^ 6 = — \E^r ^ 
where \]/ is the UOSp{l\4) Schwinger operator: 

* = {b\b 2 ,b\b\f) t . 



where the super-adjoint J is de fined in (50). It is straightforward to check that fl324) and (325) 

1 1 ■ — ' 212 J T — 

satisfy the E/OS[p(l|4) algebra (318). In the similar manner to the case of d74L we introduce 



(328) 



(329) 



Here, b a (a = 1,2,3,4) and / respectively denote four bosonic and one fermionic operators that 
satisfy 

[b Q , b^} = 6 af} , {/, /t} = 1, [b a , If] = [b a , f] = {/, /} = 0. (330) 



Square of the radius of fuzzy four-supersphere is also readily derived as 

X a X a + 2C a pe a Q p =(^\ + 3). 



(331) 



In the Schwinger formalism, the Casimir (321) is represented as 

K = XabX ab + C aP Q a e p = i (~\ (*t*)(*t$ + 3), 

a<b ^ ' 



where 



x ab = ^&r ab y. 

a 



(332) 



(333) 
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Notice that in the Schwinger operator formalism the Casimir (332) is equivalent to (331) up to 
1/2 coefficient on the r.h.s. The basis states on fuzzy four-supersphere are given by the graded 
fully symmetric representation: 



Pl, 12, 713,714) 



\mi,m 2 , m 3 ,m4) 



VnjJ 712! 713! 714! 
1 



b lT 6 2T 6 3T 6 4T , ac 



fel t^i 6 2t m2 6 3t m3 6 4t m4 yt| 



vac 



(334a) 
(334b) 



\/mil 7712! 771-3! 7714! 

where 7ti, 712, • • • , 7774 are all non- negative integers satisfying the constraint, n\ + 712 + 71.3 + 77,4 



7771 + 7772 + 7773 + 7774 + 1 
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Therefore, the dimensions of bosonic states (334e 



and 



fermionic states (334b) are respectively given by 



D B = D(n) =y{n + l)(n + 2)(n + 3), 



D F = D(n-1) 



3! 



n(n + l)(n + 2) 



and the total dimension is 



D T = D B + D F = -(n + l)(n + 2)(2n + 3). 



6 



(335a) 
(335b) 

(336) 



The bosonic and fermionic degrees of freedom (335) are respectively accounted for by the SO (5) 
symmetric basis states with the Casimir indices, n and n — 1, and hence the jV = 1 fuzzy four- 
supersphere can be regarded as a "compound" of two fuzzy four-spheres with different radii n and 
77—1: 

Sf(n)(BSf(n- 1). (337) 



S?\n) 



The eigenvalues of X§ take the following values, 

R 



d 



(n - k), 



(338) 



with k = 0, 1, 2, • • • , 2n. These eigenvalues are equal to those of X% of fuzzy two-supersphere (60) 
except for the degeneracy of the basis states at each latitude. The degeneracies at the latitude 
(338) for even k = 21 and for odd k = 21 + 1 are respectively given by 

D k=2 i(n) = (n-l + l)(l + l), (339a) 
D k=2l+1 (n) = (n-l)(l + l), (339b) 



which reproduce the total dimensions of bosonic and fermionic basis states (335) as 

n n—1 

^2 D k=2i(n) = D B , ^ D k=2l+l(n) = D F . 

1=0 1=0 



(340) 



From (331), one may find that the condition for fuzzy four-supersphere is invariant under the 
SU(4\1) rotation of the Schwinger operator ^, which is larger than the original UOSp(\\4) sym- 
metry. It may be pedagogical to demonstrate how such "hidden" SU{A\1) structure is embedded 
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in the algebra of fuzzy four-supersphere. First notice that the fuzzy four-supersphere coordinates, 
X a and @ a , do not satisfy a closed algebra by themselves, 

[x a ,x b \ = i A fx ab , [x a ,e a } = ^ a ) Pa , {e Q ,e /3 } = ^^(c 7ab ) Q/3 x afe . (341) 

a<b 



On the right-hand sides of (341), there appear "new" operators: 



X ab = ^&TabV, 9 a = ^A.*, (342) 



where T ab and D a are respectively given by (324) and 

D » = ^U). o)- <-> 

X ab and a respectively act as 50(5) antisymmetric 2-rank tensor and spinor. Commutation 
relations for these new operators can be derived as 

R R 

[X o ,0 a ] = (7a) /3a 0/3, [X a ,& a ] = ~T (7a)^a@/3, 

[X ab ,G a ] = — ( r )/ab)paQ/3, [X a b , &a] = —^{lab) /3a0 13-, 

{©a, 0/3} = y^^(C^a b ) a /3X ab , {0 a , Op] = — }XCl ab ) a pX ab , 
a<b a<b 

{G Q , P } = ^(C la ) af5 X a + ^C a pZ, (344) 
where the last equation further yields another new operator 

Z = (345) 



with 

H-( 

4 



H= U °). (346) 



Commutation relations with Z are obtained as 

[z,x a ) = [z,x ab ) = o, [z,e a ] = ~e a , [z,e a ] = ~e a . (347) 

Consequently, for the closure of the algebra of the fuzzy coordinates X a and @ a , we have intro- 
duced the fifteen new "coordinates", X ab , a and Z. In total, with the original coordinates they 
amount to twenty four operators that satisfy the SU(4\1) algebra. The basic concept of the non- 
commutative geometry is the algebraic construction of geometry, and the S'Z7(4|1) structure, the 
symmetry of the basis states on fuzzy four-supersphere, has indeed appeared as the fundamental 
algebra of the fuzzy four-supersphere. 
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7.2 U0Sp(l\4) SVBS States 



(348) 



(349) 



Similar to the U0Sp(l\2) SVBS case, the UOSp(l\4) SVBS states 46 can be constructed as 

Isvbs) = IJ(*IW ^i|4 *K r )) M l vac ) 

(id) 

= ll^bf - bfbf + 6f 6f - bfbf - r/t/t)^|vac), 

i 

where \]/(r) denotes the parameter-dependent C/05p(l|4) Schwinger operator 

and 7^i|4 signifies the UOSp{l\A) invariant matrix 

/ 1 \ 

-10 

71x14 = 1 0. (350) 

0-10 

\ -1/ 

The SO (5) spin magnitude on site i reads as 



Si = \(h?b\ + bfbj + bfbf + bfbf) = M, M-\. 



(351) 



In the following, we focus on the M = 1 UOSp{\\4) SVBS chain and its corresponding SMPS 
representation: 

L 

|SVBS) ai , afl = n(^i|4**+iW,aJvac> = (A1A2 ■ -A L )a L 



(352) 



i=i 



where A is a supermatrix given by 

/ \l,2)i V2\2,2 n 

-V2\l,l) —11,2 

A = ^n4**(r)*t(r) 



|2,3)i 
-|1,3> 



|2,4)i ^|2>i \ 
-|1,4), -^1) 



2,4) |3,4}, v^|4,4) i v^|4)i 

-^|3>i 
/ 



|M) 

—]1, 3) -|2,3)< -V2|3,3) -|3,4)< 
V-v^|l)< \^|2) 4 -V^|3) l -Vr|4)< 



with the matrix elements 



|a,a) = -^=(&°^) 2 |vac) (no sum for a), 
v2 



a,/3)^ /3 = 6 at 6 /3t |vac), 
a) = 6 a t/t h 



vac). 



(353) 



The basis states \a, a) and \a, (3) ai kp are SO (5) 10-dimensional adjoint representation while \a) are 
SO (5) 4-dimensional spinor. In total, the components of A consist of UOSp(l\4) 14- dimensional 
representation of the graded fully symmetric representation (334) for n = 2. 



46 To be precise, there exist two types of UOSp(\\4) VBS states, one of which is the tensor type (3481 and the 
other is the vector type [10]. Here, we focus on the tensor type. 
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Figure 37: The behaviors of the Schmidt coefficients and entanglement entropy (inset) of the 
5 = 1 UOSp(l\4) VBS chain. [Figure and Caption are taken from Ref. [TO].] 



7.3 Entanglement spectrum and (Z 2 x Z 2 ) 2 symmetry 

For the UOSp(l\4) SVBS infinite chain, the Schmidt coefficients are computed as 



Ab 2 = Ai 2 



Af 2 = A5 2 



A 2 



1 



A3 2 = A4 2 



8 8V25 + 16r 2 ' 



2 2V25 + 16r 2 



(354a) 
(354b) 



The bosonic Schmidt coefficients are quadratically degenerate while the fermionic one is non- 
degenerate. The behaviors of the Schmidt coefficients and the entanglement entropy, Se.e. = 
Ap log 2 Ap, are plotted in Fig. 



-4A| 



log 2 A 
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The qualitative behaviors of the entanglement 



spectra of the UOSp(l\4) VBS chain are quite similar to those of the type-I SVBS chain [see Fig. 
30 except for the quadratical degeneracy in the blue curve. 

The degeneracy of the entanglement spectra of the UOSp{\\4) SVBS chain can be understood 
based on the arguments of the symmetry protected topological order. Before proceeding to the 
case of UOSp(l\4:), we introduce the original arguments for its bosonic counterpart, the SO (5) 
VBS state |107|, 11081 1109] . Each of inversion symmetry and time reversal symmetry guarantees at 
least double degeneracy of the entanglement spectra of the 50(5) VBS states as proven by similar 
manner to the SU(2) VBS states. This is simply because each of the inversion symmetry and 
time reversal symmetry is a realization of Z 2 symmetry. The crucial difference to the SU(2) case 
is the existence of (Z 2 x Z 2 ) 2 symmetry of the 50(5) VBS states originating from the 50(5) spin 
degrees of freedom at edge [205]. The 50(5) rotation represents a rotation in five- dimensional 
space, and we "divide" the five-dimensional coordinates (1, 2, 3, 4, 5) to two three-dimensional 
subsectors, (1, 2, 5) and (3, 4, 5). The 7r-rotational symmetry around x and z axes in each sector 
generates the Z 2 x Z 2 symemtry and in total two independent sets of such discrete rotations give 
rise to (Z 2 x Z 2 ) 2 symmetry in the five-dimensional space. The group elements of (Z 2 x Z 2 ) 2 
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consist of 16 bases: 



z 2 xz 2 z 2 xz 2 



/ " v / " V 

(l,u 12 ) x (l,ni 5 ) x (l,u 34 ) x (1,1135) 



(355) 



where u a b are the 50(5) group elements of it rotation generated by the SO(5) generator a a i>: 



As the Z2 x %2 symmetry generates at least doubly degeneracy in the entanglement spectra, the 
(Z2 XZ2) 2 symmetry guarantees the four-fold degeneracy of the entanglement spectra of the SO(5) 
VBS chain. For the SUSY case, just as we have discussed in the U0Sp(l\2) SVBS case, symmetry 
transformations are attributed to those of the bosonic and fermionic sectors, 



Similar to the discussions about double degeneracy in the U0Sp(l\2) case, either of the bosonic 
and fermionic sectors brings the quadruple degeneracy to the entanglement spectra of the UOSp{\\4) 
SVBS state in the presence of (Z2 x Z2) 2 symmetry. 

8 Summary and Discussions 

We reviewed the constructions and basic properties of the fuzzy superspheres and SVBS models. 
Particularly, fuzzy superspheres and SVBS models with UOSp(l\2), UOSp(2\2) and UOSp(l\4) 
symmetries were discussed in detail. We clarified the mutual relations among the fuzzy spheres, 
QHE and VBS states were also emphasized based on the Schwinger operator formalism. It was 
illustrated that, though the SVBS states incorporate fermionic degrees of freedom, they "inherit" 
all the nice properties of the VBS: 

• Solvable parent Hamiltonian 

• Gapped bulk and gapless edge excitations 

• Generalized hidden order. 

We explicitly derived the spectra for gapped excitations (magnon and hole excitations) on ID 
SVBS chain within SMA. Physical properties of the SVBS models are qualitatively different from 
those of the other generalized VBS models based on bosonic Lie groups, in the sense that the 
SVBS states accommodate the charge sector in addition to the spin sector. In each sector, the 
SVBS states exhibit the following properties: 

• In charge sector, the SVBS states have the superconducting property (SSB). 

• In spin sector, the SVBS states show non-trivial topological order of QAFM (no SSB). 



(356) 




(357) 
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We also established the SMPS formalism for the SVBS states, which naturally incorporates the 
edge degrees of freedom to provide a powerful tool in investigating topological property. For 
practical use, the (S)MPS formalism greatly simplifies the calculations of physical quantities, 
such as the string order and entanglement spectrum. The SVBS states bear a finite string order 
regardless of the parity of the bulk superspin unlike the original VBS states. From the explicit 
calculations of the entanglement spectra of the SVBS chains, we demonstrated that there exists 
a hallmark of topological phase, the double degeneracy in SUSY entanglement spectrum. The 
degeneracy is naturally understood by invoking particular edge states in SUSY chain: Since 
SUSY relates the integer and half-integer edge-spin states, the SVBS chains always accommodate 
half-integer edge spin (as well as integer edge spin) that guarantees at least double degeneracy 
of the entanglement spectra. Consequently, the topological order is stabilized regardless of the 
parity of bulk-spin in the presence of SUSY. 

Though we focused on the SVBS states, the arguments of the present SUSY protected topo- 
logical order are applicable to general boson-fermion systems. By reformulating boson-fermion 
system, such as boson-fermion mixture cold atom system [203] with SUSY, we may apply the 
present results to discuss the stability for their topological phases. The idea of topological insu- 
lator has begun to be applied to particle theory model [204J. It may also be interesting to apply 
the present arguments to other SUSY models that are not directly related to condensed matter 
physics. 
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Appendix 

A UOSp(J\f\2K) 



Here, we summarize the basic properties of UOSp(M\2K) algebra. We denote the generators of 
the orthosymplectic group OSp(J\f\2K) as T,ab, which satisfy the following relatior 47 



■hi< 




+ 



■hi< 




Sab = 0, 



(358) 



47 The minus sign in front of in (3581 is not important for the definition of OSp(JV\2K), but added to be 
consistent with the notation of the present paper. 



87 



where the supertranspose, st, is defined in (275), lj\f denotes M x M unit matrix, and J^k 
represents the Sp(2K, C) invariant matrix 




^ = (_ lK „j- m> 

Y*AB can be expressed by a linear combination of the following matrices: 






where a and /3 are the indices of Sp{2K, C) (a, /3 = 1, 2, • • • , 2/f ) and / and m are those of 0(N) 
(l,m = 1, 2, • • • , A/"), via stand for arbitrary 2K x 2K matrices, while a a p and ai m respectively 
signify 2K x 2K and N x N matrices that satisfy 

<Hm + <?im = 0, (361a) 
cr a ^J2K + J2Kcr a /3 = 0. (361b) 

The OSp{M\2K) algebra contains the maximal bosonic subalgebra, sp(2K, C) © o(M,C), whose 
generators are <j a p and ai m . The off-diagonal block matrices T,i a are fermionic generators that 
transform as the fundamental representation under each of the transformations of Sp(2K,C) 



and 0(M, C). The o(M, C) matrices <7; m are antisymmetric matrices (361a), and then we can 
take the indices of ai m to be antisymmetric, oi m = —o m i. In the following, we consider the 
real antisymmetric matrices, the generators of o(M), with M(M — l)/2 real degrees of freedom. 



Meanwhile from the relation (361b), the generators of Sp(2K,C) a a p take the form of 

<r a p=r, V), (362) 



where k stands for an arbitrary K xK complex matrix, and similarly s and s' are KxK symmetric 
complex matrices. If the Hermiticity condition is imposed, a a p are reduced to the generators of 
USp(2K) that take the form of 

"Vtf=^ t _{*y (363) 
where h represents an arbitrary Hermitian matrix and s also signifies an arbitrary symmetric 



complex matrix. Consequently, UOSp{N\2K) generators consist of (360) whose blocks satisfy 



(361) and (363). The real independent degrees of freedom of o a p are K(2K + 1). Then, for 
usp{2K) matrices cr a /3, we can take the indices to be symmetric, a a p = ap a . Meanwhile, the real 
degrees of freedom of the fermionic generators Y>i a are 2KN. Then in total, th real degrees of 
freedom of uosp{N\2K) are given by 

dim[uosp(Af\2K)} = - (4K 2 + N 2 + 2K - N) \ 2KN = + Nf + 2K - N). (364) 
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In the present paper, instead of J 2 k (359), we adopt the following 2K x 2K matrix (which is 
unitarily equivalent to J 2 k)'- 



(365) 



n 2K 

and the following UOSp(J\f\2K) invariant matrix 



^ia 2 








o\ 





ia 2 



















V° 








ia 2 J 



n 2K 
-Iat 



(366) 



For instance, for the UOSp(l\2), UOSp(2\2), and UOSp(l\4), TLj^\ 2K are respectively given by 





( '° 


1 


o\ 




^•l|2 — 


-1 








> ^2|2 — 









-v 





-10 
0-10 
\ -1/ 



1 Kl\4 



( 1 \ 

-1 

1 

0-10 

\ -1/ 



(367) 



B Fuzzy Four-Supersphere with Higher Super symmetries 



By generalizing the construction of M = 1 fuzzy four-supersphere based on UOSp{\\£) [Sec. 7.1 
we can construct TV-SUSY fuzzy four-sphere with use of the UOSp(N\A) algebra |123j . The 
dimension of the UOSp{N\ / \) algebra is given by 



dim[uosp(AA|4)] = 10 + ^M{M - 1)|4A/" = 10 + + 7). 



(368) 



We denote the bosonic generators of uosp{M\4) as T a b = —Tb a (a,b = 1,2,3,4,5), 

Tim — -Tmi 

(l,m = 1, 2, • ■ • ,M), and fermionic generators as Fi a (a = 6\, 9 2 , #3, #4). In total, they satisfy 

[Tab,Ti a ] = (Tab)^! 1 ;^, 
[Tab, Tim] = 0, 

{r«a,r»n/3} = y^^Cja^apTgbSlm + -C^f/m, 
a<b 

{TlajTrnn] — ip/mn)plT pai 

\Tlmi T n p\ — dlnTmp ~t~ $lpTmn $mpTln ~\~ ^mn-T/p, (369) 

— 1ml (I < m) are SO(M) 
(370) 



where C is the 50(5) charge conjugation matrix (319), and ji r 
generators given by 

{.1lm)np — ^In^mp $lp$mn- 



89 



We introduce the coordinates of S, 



R 



4\2Af 
f 



18 



R 



x a = -^r a v, eg) = -*tr fa vfr, Yi 



d 



d 



(371) 



where 9 = (b±, 62, 63, 64, /1, / 2 , • • • ; /a/")* denotes the J70S7j(A/"|4) Schwinger operator whose first 
four components, b a (a = 1, 2, 3, 4), are the bosonic operators while the remaining J\f components, 



fi (I = 1,2, • • • , A/"), are the fermionic operators. The sandwiched matrices in (371) 

r a 



7a 

(V 



r 



la 



( O3+Z T a 

-(Cr a )* 
0^-/ 



lm 



o 4 

, Jlrr, 



(372) 



are the fundamental representation matrices of uosp{M\4). Here, denotes k x k zero- matrix, 



and r a are given by (326). Square of the radius of the A^-SUSY fuzzy four-supersphere is readily 
derived as 

X a X a + 2 C a pe$ef + ]T = ( ~J ) n(n + A-M), (373) 

1=1 Km=l ^ ' 

with n = ty'ty. Notice that the square of the radius of AA-SUSY fuzzy supersphere is proportional 
to n(n + 4 — M) and takes negative values for sufficiently small n that satisfy n < N — 4. For the 
positive definiteness of the radius, the SUSY number should be restricted to A/" < 4. The basis 
states on AA-SUSY fuzzy supersphere ^ are given by the graded fully symmetric representation: 



1 



-btbtbtbt 



Vh\ w. h\ k\ 1 - 1 ; 1 
1 



vac 



\m 1 ,m2,m 3 ,m i )i 1 



b\ m % m % m %rfi\^) 



\fm\\ m^. 7713! 7714! 1 
\ni,n 2 ,n 3 ,n A )i 1<i2 = = ^ =b\ ni b\ U2 b\ n3 b\ n * f} fUvac) 



1 \ 1 i3 qi h\ q2 u\ qs u\ qi A -ft t\ t\ 1 \ 

kl,g2,g3,g4)f 1 <i 2 <-<w_ 1 = ^/ gi l g2!g3!g4 | 6 l 6 2 & 3 6 4 fUUL'" /W-l l VaC > 



ki,r 2) r 3 ,r4) 



1 



where Z1+Z2+Z3+Z4 = mi+W2+™3+™4 + l = ni +722+713+714+2 



bPCCC 44/3 • • • /Jtf-i/Jtf |vac>, (374) 

• • = <7l+g2+<?3+<?4+AA-l = 



r i + r 2 + r 3 + r 4 + A" = n. Therefore with D(n) (335a), the dimension of (374) is given by 



D T = Y 



1=0 



j • D(n - I) = |(2n + 4 - AT) ( (2n + 4 - AT) 2 - I + •'W J 2 



JV-4 



(375) 



48 X a , 0q\ Yi m do not satisfy a closed algebra by themselves, and the minimally extended closed algebra including 
these operators is SU(4\Af) [T23] , 
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for n > N — 3. The degrees of freedom of the basis states (374) imply that sf 2 "^ {n) can be 
interpreted as a "compound" of lower SUSY fuzzy four-spheres with different radii: 



S, 



4\2Af 



n ) — ^2 l C m-Sp 2Af 2 \n-m) 



m=0 



Q 4\2Af-2l ( 



n) e 



I x S, 



A\2M-2l 



(n-1) 



2! 



\n-2) 



S. 



4\2N-2l 
f 

(376) 



(n 



More explicitly, 



5, 



4|W - 2 (n)0,S f 4| ^- 2 (n-l) 



^ 5 f 4 i 2Ar - 4 (n)e 



2xSp 2M ~\n-l) 



3x S, 



4\2Af- 
f 



'(n-1) 



S? 2 *- i (n-2) 



3x5, 



4|2Af-6 
f 



(n-2) 



5 



4|2A/"- 



'(n-3) 

(377) 



41 2 

For jV = 1, this reproduces the relation for 5 f (n) (61 ). 
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